APPROXIMATING TRAVELLING WAVES BY EQUILIBRIA OF 
NON LOCAL EQUATIONS 

JOSE M. ARRIETA*, MARfA LOPEZ-FERNANDEZt , AND ENRIQUE ZUAZUA * 

Abstract. We consider an evolution equation of parabolic type in R having a travelling wave 
solution. We perform an appropriate change of variables which transforms the equation into a non 
local evolution one having a travelling wave solution with zero speed of propagation with exactly 
the same profile as the original one. We analyze the relation of the new equation with the original 
one in the entire real line. We also analyze the behavior of the non local problem in a bounded 
^— H interval with appropriate boundary conditions and show that it has a unique stationary solution 

T— H which is asymptotically stable for large enough intervals and that converges to the travelling wave 

as the interval approaches the entire real line. This procedure allows to compute simultaneously 
the travelling wave profile and its propagation speed avoiding moving meshes, as we illustrate with 
several numerical examples. 

Key words, travelling waves, reaction-diffusion equations, implicit coordinate-change, non- 
local equation, asymptotic stability, numerical approximation. 
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p ^ 1. Introduction. We address the problem of the analysis and effective compu- 

tation of travelling wave solutions emerging from parabolic semilinear equations on 
the real line: 

(-H 

f Ut{x,t) = Uj:x{x,t) + f{u{x,t)), — CXD < X < -l-CXD, < > 0, 

< (1.1) 
g [ u{x,0) = Unix). 

We assume that f G with /(O) = /(I) = 0, so that u = and u = 1 are stationary 



I solutions of (1.1). Under these assumptions, if the initial data uq is piecewise contin- 

^ uous and < ^^o < 1, there exists a unique bounded classical solution u{x,t) defined 

for all i > and, due to the maximum principle, < u{x, t) < \ for all x, t. 
*^ A travelling wave is a solution of the type u{x, t) = $(a; — ct) where the function 

$ is the profile of the travelling wave and c is the speed of propagation of the wave. 
For instance, if c > (resp. c < 0) the solution will consist of the profile x — $(a;) 
'- 1 travelling in space to the right (resp. left) with speed |c|. Proofs of the existence of 

this kind of solutions can be found in [TJ [13 [H] among others. 

The asymptotic profile $, when it exists, will have finite limits at ±oo, either 
$(—00) = 0, $(00) = 1 or $(—00) = 1, <i>(oo) = 0. In the first case $ will be a 
. solution to 

>< r $"(o + c<i>'(e) + /(<i>(0) = 0, -oo<c<+oo, 

H { (1-2) 

C3 1 < $ < 1, $(-00) = 0, $(+00) = 1, > 0, 
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that is called a [0, Ij-wave front. The monotonicity condition on the profile $ is not 
a restriction but rather an intrinsic property of the travelling wave profiles, as it is 
shown in [T31 Lemma 2.1.]. 

Note also that, by simply making the change of coordinates £, -t> — ^, to every 
pair ($,c) with $ a monotone increasing [0,l\-wave front corresponds a monotone 
decreasing [1, 0\-wave front with propagation speed — c. Thus, all what follows applies 
to monotone decreasing solutions as well. 

These profiles are well known to have the property of attracting, as i — >■ oo, 
the dynamics of a significant class of solutions of the Cauchy problem (jOJ, see for 



instance Theorem 2.1 below (from [13]). In [T51 it is proven that if / is of bistable 
type satisfying 

/(0) = /(l)=0, 

/'(O) <0,/'(l) <0, (1.3) 
3a e (0, 1), s.t /(u) < 0, for ue(0, a), f{u)>0, for ue{a,l). 



then, for a certain set of initial data uq, the solution u{x,t) of (1.1 1 evolves into a 
travelling wave ^{x — xq — ct), for a certain € K depending on the initial datum 
Uq, i.e., 

\u{x,t) - ^{x - xo ~ ct)\ —i' 0, as t->oo. (1.4) 



Further, the convergence in (1.4) is shown to be uniform in x and exponentially fast 
in t. 

From the numerical point of view, one of the main difficulties in the approximation 
of these asymptotic solutions $ and their propagation speed c is the need of setting 
a finite computational domain. While the solution u evolves into $, it also moves left 
or right at velocity c and it eventually leaves the chosen finite computational domain. 
A natural approach is then to perform the change of variables u{x,t) = v{x — ct,t), 
so that the resulting initial value problem for v converges to a stationary solution, i.e, 
f ( — >■ as i — ?> oo. However, in general, the value of c is not known a priori. 

The issue of having a priori characterizations of the velocity of propagation c 
then plays an important role from a computational viewpoint. It has been addressed 
in a number of articles. For instance, in |28j explicit mini-max representations for 
the speed of propagation c are provided. Unfortunately, the expressions in ^S] are 
difficult to handle in practice for the effective computation of c and $. 

In the present paper, we follow and further develop the approach introduced in 
[5], where a new unknown 7(<) is added to the problem, to perform the change of 
variables 



u{x,t) ^ v{x - ■y{t),t). (1.5) 

Then, v satisfies the equation 

Vt{x, t) = v,^^{x, t) + ^'{t)v^{x, t) + f{v). 

Our goal is to determine 7 = "f{t) a priori so as to ensure that, as time evolves, 
it converges to the asymptotic speed c of the travelling wave. Of course, in order 
to compensate for the additional unknown 7, one has to add a so called "phase 
condition", that is, an additional equation linking v and A := 7'. Two different 
possibilities were proposed in [5]. One of them consists on minimizing the L^-distance 
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of the solution w to a given template function v{x), which must satisfy u — $ e H^{M.). 
This approach leads to a Partial Differential Algebraic Equation of the form 



+ Xvx + f{v), 



'CO < X < oo, t > 0, 



= {v' ,v — v), 
v{x,0) = uo{x), 



(1.6) 



where (•,•) denotes the inner product in L^(R). This system is locally (in time) 
equivalent to the original one (1.1 1, provided that the implicit function theorem can 
be applied to the equation 

$(7,i) {v',u{- + j{t))~v)^0, 

to obtain ^{t) and v{x,t) — u{x + 'y{t),t), once u is given, see [3 Theorem 2.9]. The 
approximation properties of (1.6) and its numerical discretizations have been studied 
in detail in [551 [23 [27] . This approach is useful as long as a suitable template mapping 
V, close enough to is available and the initial data uq belongs to $ + H^{R), too. 



From the computational point of view, the semidiscretization in space of ( 1.6 ) leads to 



a Partial Differential Algebraic Equation of index 2. This means that it is necessary 
to differentiate twice the algebraic constraint in order to eliminate it and get the 
underlying Partial Differential Equation, see for instance [T71 Chap. VII]. 

A more global approach, also proposed in W, is obtained by minimizing ||wt||2- 
This yields the augmented system 



+ Xvx + f{v), 



-oo < X < oo, t > 0, 



= \{Vx,Vx) + {f{v),Vx 

v{x,0) = uo{x), 



(1.7) 



which is a Partial Differential Algebraic Equation of index 1. Assuming that the 
products above are well defined and ||fi:(-, t)||2 ^ 0, t > 0, the phase condition 



X{Vx,Va:) + {f{v),Vx) 



yields directly 



lit) 



{f{u{-,s)),Ux{-,s)) 



ds 



' {f{v{;s)),vA;s)) 



lo {Ux{-,s),Ua:{-,s)) Jq (w^; (' , s) , (• , s)) 

and one gets the following nonlocal semilinear equation 

{f{v),Vx 



Vt = Vxx 

{Vx,Vx) 

v{x, 0) — uq{x). 



Vx + f{v), — oo < X < oo, t > 0, 



(1.9) 



Definition (1.8) is in fact quite natural if we notice that, after multiplicating by $' in 

(1.10) 



(1.2) and integrating along the real line, we obtain 



($',$') 
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Thus, the equation for $ can be written in the form 



(1.11) 



Equation (1.9) is a time evolving version of ( |1.11[ ). It is natural to expect that, as 
t oo, system (1.91 wih yield both the speed of propagation and the profile of the 
travelling wave. 



We notice that (1.10) is equivalent to 



F{1) 



(1.12) 



with 



Fiu)^- / fis)ds. 
Jo 



(1.13) 

This observation leads in a natural way, to the following alternative definition of j{t) 



in (1.5) 



Fil) 



lit) = - / , , ' ds, 
which yields the nonlocal semilinear problem 

m 



(1.14) 



{Vx,Vx) 

v{x, 0) = uq{x). 



Vx + fiv), — oo < X < oo, t > 0, 



(1.15) 



We will show that ( 1.15[ ) enjoys in fact similar properties to those of (1.9) and, to our 
knowledge, it has never been used in practice to approximate travelling waves and 
their propagation speed. 

In [S], the use of (1.8) or, in other words, (1.9) or (1.15), is regarded to be 



particularly useful near relative equilibria and good numerical results are reported. 
Moreover, observe that analyzing these two equations requires less a priori knowledge 
of the asymptotic state than analyzing (1.6) (no template function in $ + iJ^(K) is 
required), and it allows to consider more general initial data. Moreover, equation 
(1.9) is simpler to approximate numerically than (1.6). However, to our knowledge. 



no rigorous asymptotic analysis of the modified equation (1.9), (1.15) seems to be 
available. In fact, in later works by the authors [51 |251 [37], where the effects of 
discretizations are also taken into account, this kind of phase condition is not analyzed. 



focusing only on the extended system (1.6). 

In this work, we set necessary conditions for the well-posedness of these two 



new initial value problems (1.9), (1.15) and analyze the relation between these two 



problems and the original one (1.1). We will prove that both problems (1.9), (1.15) 



have the one parameter family of travelling waves with c — speed of propagation 
^{x — a), a ^ M, ("standing wave") with the same profile $ of the original equation 
( |l.l[ ). Moreover, under appropriate assumptions on the initial data uq, we prove that 
j'{t) — c, as i — oo (and we recover the speed of propagation of the travelling wave 
of the original problem), and the solutions of (1.9), (1.15) converge exponentially fast 
to one of these standing waves. 



Once the modified problems (1.9), (1.15) are understood and shown to converge 



to an equihbrium state with the same profile $ as the travelling wave for (1.1), the 



problem of its numerical approximation arises naturally, but can be addressed quite 
more easily because, now, one does not need to address the issue of moving the frames 
as time evolves. To this end it is necessary to truncate the spatial domain and add 
some reasonable artificial boundary conditions. This motivates the analysis of prob- 



lems ( 1.9 ), ( 1.15 ) in a bounded spatial interval (a, b) with certain "artificial" boundary 
conditions. We have chosen non homogeneous boundary conditions of Dirichlet type 
which emulate the behavior of the travelling wave in the complete real line, that is. 



^^(1) 



''"\L^(a,b) 

v{a,t) = 0; v{b,t) = 1, t > 0, 
v(x,0) — uq{x), X G [a,b]. 



+ fiv), X e (a, 6), t > 0, 



(1.16) 



Observe that when restricting both equations (1.9), (1.15) to a bounded interval and 
impossing v{a) = 0, v{b) = 1 we obtain in both cases the very same equation, which 
is the one given above in (1.16). 



We analyze equation (1.16) and show that with the nonlinearity / satisfying 
(1.3) we have a unique stationary state $(a,b) with < $(a,6) ^ 1- Moreover, this 
stationary state, when normalized so that $(a.6)(0) = 1/2, will converge to the profile 
of the travelling wave of equation ( 1.9 ), ( 1.15 ) as (a, b) (— oo, +oo) (see the details in 
Section 4.3 ). We also analyze the stability properties of this stationary state. In order 



to accomplish this, we will need to analyze the spectral properties of the linearization 



of (1.16) around the stationary state, which means to analyze the spectra of the 



"nonlocal operator" 



Lw 



F{1) 



1$; 



{a,b)\\L^ia-b) 

2F{1) 



Wx + f'{^ia,b)i^))w 
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\L^(a,b) 



(a,fc) / W{x)<^[^,,^{x)dx. 



This task is not a simple one. There are in the literature several works which 
analyze the spectra of operators of the type above, see HOI HH HH US] j but none of 
them are conclusive enough to characterize it completely in our case. 

Nevertheless, we will be able to show that cr{L) C {z £ C,Re(z) < —K{a,b)} 
for certain n^a, b) > when the length of the interval is large enough (that is, for 
6 — a — +00). We will obtain this result via a perturbative argument, viewing the 
operator i(a,b) a perturbation of the operator on (a, b) = M. that is 



Lw = Wx 



-Wx^ f{<^{x))w 



$' / w'{x)'^'{x)dx, 



2f(l) 

for which the spectra is easier to characterize, since $' is the eigenfunction associated 
to the eigenvalue or the operator Loo- 

We will conclude in this way the asymptotic stability of the stationary solution 
^(a,6) for large enough intervals (a, 6). 

The analysis in the present work is performed for non-homogeneous Dirichlet 
boundary conditions. This choice is justified since, somehow, it imitates the behavior 
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of the traveUing wave for large enough intervals. Nevertheless, other boundary con- 
ditions may be suitable to approximate the travelling wave although the dynamics 



of system (1.16) with these other boundary conditions may differ from from the case 



treated in this paper. As a matter of fact an analysis of the differences and similarities 
for different boundary conditions will be important and will be carried out in a future 
work. 

Let us notice that the idea of performing a change of coordinates so that the 
front of the asymptotic profile $ remains eventually fixed in space appears also in 
[21], where a different change of variables is considered. This alternative change of 
variables is based on the fact that it also holds 



(1.17) 



Formula (1.17) is readily obtained after integration along the real line in (1.2) and 
leads to the change of coordinates u{x,t) = v{x — 7(t)) with 



lit) 



f{u{x, s)) dx ds. 



The convergence of the solutions of the resulting equation to an equilibrium is not 
proved in [21]. In this paper, we do not study this particular change of variables 
although we expect that by adjusting the techniques we develop here we will be able 
to obtain similar results. As a mattter of fact, we regard the analysis included in this 
paper as a general technique that, with the appropriate adjusments to the different 
possible changes of variables which make the velocity of propagation implicit in the 
equation, will yield in an effective way both, the speed of propagation and the profile 
of the travelling wave. 

The paper is organized as follows. Sections [2] and |3] are devoted to the deduction 
and analysis of the problem in the whole real line. In Section [2] besides recalling the 
result on existence of travelling waves from |13j . we also obtain several important 



estimates of the solutions of the original problem ( 1.1 ) when the initial condition uq 
satisfies Mq G and d^u^ € Lp(R), 1 < p < oo. 

The change of variables that leads to the modified problema (|1.9|), ( |1.15 ) is con 



sidered in Section [3] where we establish a fundamental relation between these new 
problems and the original one ( 1.1 ). The estimates obtained in Section [2] are used in a 
crucial way in this section. We show the asymptotic stability, with asymptotic phase, 
of the family of travelling wave solutions of the nonlocal problems ( |1.9[ ), (1.15), see 
Theorem 13.31 



The next two sections, Section [4] and Section [5] are devoted to the nonlocal prob- 
lem in a bounded interval. In Section [4] we obtain the existence and uniqueness of an 
stationary solution of problem ( 1.16[ ) and show that the stationary solution converges 
to the profile of the travelling wave solution in the entire real line. In order to accom- 
plish this, we will need to perform a careful analysis of the behavior of the associated 
local problems in a bounded domain, see Subsection |4.1[ and then to relate the results 
obtained for the local and nonlocal problem, see Subsection 4.2 The convergence of 



the stationary states to the travelling wave is obtained in Subsection 4.3 In Section 



[5] we show the asymptotic stability of the stationary state of the non local problem 
in a bounded interval. We analyze first the properties of the spectra of the linearized 
non local equation in a bounded interval, see Subsection |5.1[ and in the complete 
real line, see Subsection 5.2 In Subsection 5.3 we obtain the asymptotic stability of 
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the stationary states for large enough intervals. This is obtained through a spectral 
perturbation argument. 

Finally in Section [6] we include several mumerical examples which illustrate the 
efficiency of the methods developed in this article to capture the asymptotic travelling 
wave profile and its velocity of propagation. 

Ackov^rledgements. We would like to thank Erik Van Vleck, Pedro Freitas, 
Wolf-Jiirgen Beyn and Nicholas Alikakos for the different discussions we had with 
them and for pointing out several important aspects of the subject of this paper. 



2. Estimates for the original problem (1.1 ). We start reviewing some of the 



results in [13^ and [19] on the existence and behavior of travelling wave solutions of 



(|1^. 

The following theorem of |13| ensures the existence and uniqueness of an asymp- 



totic travelling front for the original problem (1.1) under quite general assumptions 
on the initial data uq. 



Theorem 2.1. Let f G C^[0, 1] satisfying (1.3). Then there exists a unique 



(except for translations) monotone travelling front with range [0,1], i.e., there exists 



a unique c* and a unique (except for translations) monotone solution $ of (1.2). 
Suppose that uq is piecewise continuous, < Uo{x) < 1 for all x G M, and 

liminf ito(a;) > a, limsupuo(2:) < a- (2-1) 



Then there exist Xq G M, K,uj > 0, such that the solution u to (1.1 1 satisfies 

\u{x,t)-^{x-c*t-xo)\<Ke-'^\ xgM, t>Q. (2.2) 

Furthermore, c* > (resp. < 0) when F{1) — — /(s) ds > {resp. < 0). 

Summary of the Proof of Theorem |2.1| in [13] . For c* in the statement of 
Theorem [211 set 



w{x,t) ^u{x + c*t,t), (2.3) 

which fulfils 

{Wt{x,t) = Wxx{x,t) + C*Wx{x,t) + f{w{x,t)), — CX3 < x < oo, t > 0, 
(2.4) 
w{x,0) = uo{x). 

The proof is then based on the construction of a Liapunov functional for equation 



(2.4). The main tools are a priori estimates and comparison principles for parabolic 
equations ]i6j Theorem 4 of Chapter 7 and Theorem 5 of Chapter 3] . The following 
two Lemmas are important intermediate steps in this construction and will be used 
in Section [31 

Lemma 2.2. Under the assumptions of Theorem \2.1\ there exist constants Xi, 



X2, qo and pL, with qa, p, > 0, such that 

$(2: - a;i) - goe"'"* < w{x, t) < <I>(x - X2) + goe"'"*. (2.5) 
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The following lemma provides asymptotic estimates for the derivatives of w. 
Lemma 2.3. Under the assumptions of Theorem \2. 1\ there exist positive constants 
a^fi and C with a > |c*|/2, such that 



\l~w{x,t)\, \Wx{x,t)\, \Wxx{x,t)\, \wt{x,t)\ 

< C(e-('=*/2+<T)x^g-Mt)^ a;>0, i>0; (2.6) 
\w{x,t)\, \w^{x,t)\, \wxx{x,t)\, \wt{x,t)\ 

< Cie^"-"'''^^'-' + e->''), a;<0, t>0. (2.7) 



Remark 1. Although the result stated in Theorem\2.1\ is very general in terms 



of the initial data in (1.1 1, its proof yields little hint about the exponent u) in the 
exponential estimate (2.2 1. In this sense the study accomplished in \19jl is clearer. 
Following a different approach, the convergence result (2.2 1 is also proven in \1S^ . 



although for a less general class of initial data. Once an equilibrium $ for (2.4 1 is 



shown to exist, the uniform convergence of w to a shift of $ is obtained by analyzing 



the linearization about $ of the equation in (2.4). More precisely, the spectrum of the 
operator 



Lw :— w" + c*w' + /'($)w, 



oo < X < oo. 



(2.8) 



is considered. By JlfA Theorem A. 2 of Chapter 5], the essential spectrum of L lies in 
Re z < —13 with 



/3 = min{-/'(0),-/'(l)}>0. 



(2.9) 



The rest ofa{L), i.e., the set of isolated eigenvalues of L of finite multiplicity, is also 
shown to be negative but the eigenvalue 0, which turns out to be simple. This, by \1!A 
Exercise 6 of Section 5.1], yields the exponential rate of convergence in (2.2). The 
rate of convergence can be taken as any uj < ujq where 



ujo = min{/3,7}, 



(2.10) 



where —7 < is the spectral abscissa, i. e. the largest real part of any non zero 
eigenvalue of L. In fact, this analysis of cr(L) yields the asymptotic stability with 
asymptotic phase of the family of equilibria 

{^{■-xo):xoeR} 



of (2.4) 



We finally notice that from the proof of Lemma 2. 2 accomplished in fT^ it follows 

and Lemma\2.3\ can be chosen as close to /3 as we 



that the constant fj, in Lemma 2.2 



wish. This implies that we can choose any fi satisfying 
The above bound will be used in Section\^ 



(2.11) 



We next show an existence and uniqueness result for the original Cauchy problem 



( 1.1 ) in the spaces 



W^i'f (M) = {u e Wlim : dxu G LP{R)}, 



(2.12) 



for 1 < p < oo. This result is slightly more general than what we strictly need to 

0. Letl<p<oo 



ensure the well-posedness of (1.9) 



Proposition 2.4. Let f e Ci(]R,]R) satisfying /(O) = /(I) 
and uq e L°°(M) f] W^i^P(M) with < uq < 1 a.e. x€R. Then, 

(i) there exists a unique mild solution u E L°°([0,oo) x M) P| C((0, oo); Vl^^'P(M)) of 
the Cauchy problem (1.1 1. Moreover, this solution satisfies < < 1; is a classical 



solution for t > and has the following regularity u £ C(0, go, C^'''(M)), for all rj < 1. 
(ii) In case p = 1 and when the function f satisfies [J. jji and the initial condition uq 



satisfies (2.1), then there exists C > such that \\ux{- ,t)\\i < C, for all t > Q. 



(Hi) In case p = 2 and when the function f satisfies {1.3) and the initial condition Uq 
satisfies (2.1), then there exists a (5 > such that ^)||2 > /? for all t > 0. 



Proof, (i) For initial data uq € L°°(M) the existence and uniqueness of mild 
solutions in L°°((0, oo), L°°(K)) holds by the variation of constants formula below 
and standard fixed point arguments, 

u{t) = G{-,t)*uo+ [ G{t~ s)* f{u{-,s))ds, t>0, (2.13) 
Jo 

where G — G{x,t) — (int)^^/^ exp{—\x\'^ /At) is the heat kernel and * denotes the 
convolution in the space variable. Then, by the regularization properties of the heat 
kernel, this solution belongs to G{{0, oo),W''°° (R)), for ah s < 2, see The 
embedding M^*'°°(M) ^ C^'''(E) for 77 < s — 1 implies the regularity result. 
On the other hand let us consider the initial value problem 



qt = Qxx + f'{u{x,t))q, 
q{x,Q)=dxUG{x) e LP{1 



-00 < a; < 00, t > 0, 



formally solved by the space derivative of u. Since q{x,0) e LP{1 
unique solution q e C([0, oo),LP( 



(2.14) 



we have a 



of (2.14). In fact, q is given by (2.131 with 



f'{u{-, s))q{-, s) instead of f{u{-, s)). The well-known L'' estimates for the heat 

equation in K^, namely 



\\G(;t)*^\\,<Ct- 



1 < p < q < 00, 



imply that also q e C(0, 00, i°°(]R)). But then it is 



q{x,t) ^ Ux{x,t) ^ G{-,t) * da:UQ^ / G{t- s)*{f'{u{-,s))ux(,s))ds, 



(2.15) 



t > 0, 
(2.16) 



by the uniqueness of solutions of (2.14) in C(0, 00, L 



Since /(O) = and /(I) = 0, both functions u = and u = 1 are strong solutions 



of (1.1 ). Using standard comparison arguments, we have that if < mq < 1 then any 



possible solution starting at uq will lie between this two constants functions. 

(m) From ( |2.2[ ) and i) we have that the solution u to ( |1.1[ ) approaches a travelling 
wave solution and Ux{-, t) e L^(]R) for every t. We next show the uniform boundedness 



in time of |j'Ua;(-, i)||i. This is equivalent to bound ||wa;(-, i)||i, for w in (2.3). To this 
end, we consider h = Wx, which satisfies the equation 



hxx + c*hx + f'{w)h, 



~oo < a; < 00. 



(2.17) 



Applying that /' is continuous, that both w, $ e [0, 1], and (2.2 1, we can estimate 

\f'{w{x,t))~f'i<i>{x~xo))\<C\w{x,t)~<i>{x-xa)\<Ke-'^*, xeR, t>0. 

This, together with the hypotheses /'(O), /'(I) < 0, imply the existence of L > and 
to > large enough, and /3 > 0, so that 

f'{w{x,t)) < f'{<^>{x-xo))+Ke-'^* <-l3 <0, for |x| > L, t > to- 



Multiplying formally in (2.17) by the sign of h, sgn(/i), and integrating in {x G 
> L} gives, for every t > Iq, 



d_ 
di 



\x\>L 



\h{x,t)\dx= / {K^sgn {h) + c*\h\^)dx+ / f {w{x,t))\h{x,t)\ dx 



\x\>L 



|2;|>L 



This yields, by Kato's inequality (see |20l ) and applying estimates in Lemma 2.3 
to h — Wt and h 



d 
di 

< 



''X ^XXl 



\h{x, t) \ dx 



\h\^^dx+\c*\{\imsnp\h{A'I,t)\ + \h{±L,t)\)~ (3 / \h{x,t)\dx 



\x\>L 



A/-j-±oo 



\x\>L 



< limsup{\K{M,t)\ + \h.A±L,t)\) + C - /3 / \h{x,t)\dx 



M— >-±oo 



\x\>L 



< -13 \h{x,t)\dx + C, t > to, 

J\x\>L 

where the constant C is independent of t and L. Thus, setting 



9{t) 



\h{x, t)\ dx, 



\x\>L 



multiplying the above inequality by e^* and integrating from to to t, we obtain 

git) < e-^(*-'«)5(io) + ^(1 - e'^^*-*")) < A, t> to. 
Finally, we apply again Lemma 2.3 to estimate 

j \h{x,t)\dx < A+ i \h{x,t)\dx < A + 2L sup \h{x,t)\<C, for aU < > to- 
Jvl J\x\<L xe\-LM 



For t e (0,to] we can bound directly t)|| i in (2.16) and apply Gronwall's ine- 

quality. 

The above argument can be formalized by multiplying in (2.17) by h\h\P^'^ with 
p > \. In this way we can get an estimate for \\h\\p which turns out to be independent 
of p and then take the limit as p — )• 1. Another possibility is to consider a Lipschitz 
regularization of sgn(ft,). 

(iii) From (2.2) we have that the solution approaches a travelling wave solution 
and therefore, liminft_j._|_oo ||u2:(t, •)||2 > 0, which implies that there exists a Ti and 
Pi with \\ux{-,t)\\2 > /3i for aU t > Ti. 
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On the other hand, if there exists some time < T < Ti such that ||wk(-, T)||2 — 
then u{-,T) is a constant function and therefore u{-,t) is a constant function for all 
t > T. To see this we just use the uniqueness of solutions and the fact that if the 
initial condition is a constant function, then the solution is a constant function in 
space for all forward times. Hence, i)||2 > for all t G [0,Ti] and since this is 

a compact interval and the function t — )■ i)||2 is continuous, then there exists a 

/32 > such that i)||2 > ^2 for all t £ [0, Ti]. This shows the last part of the 

proposition. □ 

Remark 2. Assuming further that /' is Lipschitz continuous, it is possible to 



prove i) of Proposition 2.4 by using standard fixed point arguments in the space 
L°°(K) n W^i'P(K). 



Remark 3. Under the assumptions of Theorem 



tion 



2.4 



2.1 



and dxUQ e 

the solution u of (|1.1|) approaches a solution $ of (|1.2|), as i — > oo, and, by Proposi- 



$ must fulfil besides the added integrability condition $' e 
We show below that this is consistent with the following properties of such a travelling 
front <i>: 

(i) $' e iy^(M). This is a direct consequence of the fact that $' > 0, so that 
|$'| = and $(±00) < 00. 

(ii) $' e L^{M 

obtaining 

jb jr jk 



To see this, we multiply in (1.2) by and integrate along 



1 r d 

F($(oo)) - F(<I>(-oo)) = F{\) < OD, 



where F is defined in (1.131 and we used that, necessarily, $'(±cx)) — (see 

m)- 



(iii) We also notice that $" G L'^(M), too. This follows again from (1.2), now after 
multiplication by which leads to 



+00 



/+00 
/'(<i>)($')'rfe<oo, 
'OO 



since, by hypothesis, / e C^, $ is bounded and we have shown that $' G 
3. The nonlocal problem in the entire real line. We now turn to the non 



local Cauchy problems (1.9) and (1.15) 



The following result shows how the well-posedness of equation (1.9) depends on 



the properties of the solution u to the original problem ( 1.1 ). 



Proposition 3.1. Assume that the initial data uq in (1.1) is piecewise continuous 



and < mq < 1. Let u be the unique classical solution to (1.1). Assume further that 

(i) u^{-,t) e Li(M)ni^(K); for every t > 0. 

(ii) there exists (3 > such that ||ux(',^)||2 > /3 for all t >0. 
Then 



v{x,t) := u{x + ju{t),t), xeR, t>0 
11 



(3.1) 



with 



lu{t) :-- 



{ua:{-,s),u^{-,s)) 



ds. 



t > 



(3.2) 



is well defined and is a classical solution of ( |1.9[ ). 

Proof. Due to assumption (m), u{-,t) is not constant in space for all i > and 



_ {f{u{;t)),U.^{;t)) 

{u^(-,t),u^{-,t)} 

defines a bounded and continuous mapping of t. The integral in the scalar product 



of the numerator in (3.3) is convergent, since u e L°°(R), then f{u) € L°°{R), and, 
by hypothesis (i), € L^(]R). By hypotheses {i) and (iz), the scalar product in the 
denominator in (3.3) is also finite and strictly positive. The continuity follows from 



the fact that m is a classical solution. Thus, hi (3.2) is well defined and so is v in 



(3.1) 



From the invariance with respect to translations of the integral in the whole real 
line, we easily get that for every t>Q fixed, {f{v{-,t)),Vx{-,t)) = {f{u{-,t)),Ux{-,t)). 
In a similar way, |jwa;(-, = \\ux{-,t)\\l. Thus, A„(i) = A,„(<), 7„(t) = jy{t), 



v{x,t) = u{x + ^^{t),t), 



and clearly v fulfils (1.9) 



(3.4) 



□ 



The analysis of (1.151 follows the same steps as the analysis of (1.9) and is, in 



fact, simpler. Thus, we only state here the corresponding result for (1.151 



Proposition 3.2. Assume that the initial data uq in ( 1.1 1 is piecewise continuous 



and < uq < 1- Let u he the unique classical solution to (1.1). Assume further that 

(i) u,i-,t) e l2( 



), for every t > 0. 
(a) There exists j3 > such that \\ux{-,t)\\2 > (3 for all t >0. 
Then, 



v{x,t) := u{x + ju(t),t), 



with 



lu{t) 



{u,^{-,s),u,^{-,s)) 



X e 



ds, 



t > 



t > 0, 



(3.5) 



(3.6) 



and F in (1.13) is well defined and is a classical solution of (1.15 I. 



( 1.15 ), respectively, can be reduced to a further study of the original Cauchy problem 



Propositions 3.1 and 3.2 imply that the study of the well-posedness of (1.9) and 



( 1.1 ). In fact, all we need is to ensure that the solution u of ( 1.1 1 fulfils assumptions (?) 



and (m) of Proposition 3.1 or Proposition 3.2 As summarized below, this is provided 
by Proposition |2.4[ 



We are now in the position to prove the main result of this section and one of the 
main results of this paper. 

Theorem 3.3. Under the hypotheses of Theorem 2.1 and assuming further that 
Uq G W^'^(M.)f]W^''^{M.), the augmented problem (1.9) is well-posed and its solution 
V is given by (3.1 )-( [3?2| . 

Let any u) < ujq, wq being defined as in (2.10). Then, there also exist x* € 
positive constants Ci, C2, such that 
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and 



(i) for c* the propagation speed in Theorem 2.1 and in (3.3) it holds 

|A„(t) -c*| < Cie-"^*, t>Q. 



(3.7) 



(a) For ^ the unique (except for translations) solution to (1.2), we can estimate 
\v{x,t)-<^{x-x*)\<C2e-'^\ xe. 



t > 0. 



(3.8) 



The above result validates the change of variables (3.1)-(3.2), shows that Ai, in 
(3.3) converges to the asymptotic speed c at an exponential rate and provides the 
analogue to Theorem 2.1 for v, since (3.8) is equivalent to 

\u{x,t)~<P{x--f{t)-x*)\ <C2e-"* 



X e 



t > 0. 



(3.9) 



Moreover, the rate of the exponential convergence is the same as the one derived in 



Proof. By applying Proposition 2.4 with p = 1 and p = 2, we obtain that 
Ua:i-, t) g L^{M) n for u the solution to (O). Then, by {Hi) of Proposition 



2.4 



the assumptions of Proposition 3.1 arc fulfilled and the well-posedness of (1.9) follows 
in a straightforward way. 

In order to prove the convergence results (i) and (m), we also need some integra- 
bility properties of u^x- More precisely, we will use that 



u,j:x{-,t) e L^{M) and \\uxx{-,t)\\LMs.) < C, for alH > 1. 



(3.10) 



The estimates in (3.10) can be proved by considering the new variable q = u^, which 
satisfies the initial value problem (2.14). By the regularization properties of the 
equation in (2.14), see for instance 19 , it is possible to show a bound of the type 
hx{-,t + 1)||li(r) < C||g(-,t)||ii(K). This implies ( |3.10[ ), since by (m) of Proposi- 
tion 2.4 ||q(-,t)||i is uniformly bounded in t. 

(i) By using that the inner product in L^(M) is invariant under translations in 
the space variable and (Hi) of Proposition 2.4 we can estimate 

1 



|A.(i)-c*| 



{f{v),v,)+c*\\v,\\i <C {f{v),v,)+c*\\v,\\i 



for some C > 0. Then, for w defined in (2.3), $ = ^{x — xo), with a;o in (2.2), and 
formula (1.10), it follows 

{f{v),v^) +c*||-i;^||2 = {f{w),w^) +c*||w^||2 

= {f{w),w^)+c*\\w,\\l - (/($),$') 11^ 



On one hand, using that /' is continuous, w, $ e (0, 1), and by (m) of Proposition 2.4 
it is ||u;a;(-, t)||i = ||u2;(-, i)||i < C, for all i > 0, we can estimate 



\{f{w)~ fm,w,)\<C\\w-^oo\\w^\\l<Ce- 



^OT uj in (2.2). 
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Moreover, with estimates in Lemma |2.2[ Lemma |2.3| and the fact that $ ap- 
proaches its hmits at ±00 exponentially fast, there exists a > such that 



L^oo 



< limsup 

L— ^00 



f{^){w.,-<i>')dx 



fmw{;t) - $)|^_^ + / /'(<i>)<i>'(^« - $) 



< lim |C(e-'"^ + e'^*)| + !(/'($)$', 

< Ce-"* + C\\w - $||oo < Ce-"' + Ce"^\ 



In a similar way. 



— fj,t 



(3.11) 



where we used that w^x and $" are in _L-'^(]R), that ||u>j;2^||i — \\uxx\\i and (3.10). 
This implies that 

\Xv(t)-c* \ < Ce-"' + Ce-'^' 

Observe now that from Remark [T] we can choose both ^,0; < wq but arbitrarily 
close to ujQ (defined in ( 2.10[ )). Hence if uj* < ujq, we may choose uj* < /i,a; < wq 
which implies Ce~'^* + Ce^ < Cie~"'* and this shows ^jl. 

(ii) The estimate (3.71 implies the convergence of the integral 

/"OO 

/ {Xv{s) - c*)ds ^ A, for some A g M, 
Jo 

since it is absolutely convergent. Moreover, 



|7„(<)-c*t-A| = 



{\v{s) — c*) ds 



■ds^Ce 



-U}*t 



Then, setting x* = xq — A, for xq in (2.2 1, and applying Lemma 2.3 Theorem 2.1 
and (2.111 it follows 

\v{x,t) - <P{x - x*)\ = \w(x + A + 7t,(t) - c*t - A,t)~ ^(x + ^ - a;o)| 

< \'w{x + A + 7^,(t) - c*t - A,t)- w{x + A,t)\ + \w{x + A,t) - ^{x + A - xo)\ 

< sup(|w:r (a;, t)\) |7„(i) - ct - A\ + \w{x + A,t) - <^>{x + A - xo)\ 

which concludes the proof of the theorem. □ 



The analogous result holds also for problem (1.15). 



Theorem 3.4. Under the hypotheses of Theorem \2.1\ and assuming further that 
uq g T4^"'^'^(IR), the augmented problem (1.15) is well-posed and its solution v is given 
by 
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Let any lo < ujq, wq being defined as in (2.101. Then, there also exist x* € M and 
positive constants Ci, C2, such that 



i) for c* the propagation speed in (1.2) and Aj,(i) — ^y{t) for "fv{t) in (3.6| it holds 

\Xv{t)-c*\<Cie-'^*, t>0. 



(3.12) 



ii) For ^ the unique (except for translations) solution to (1.2 I, we can estimate 



\v{x,t) - <S>{x - x*)\ <C2e~ 



X e 



t > 0. 



(3.13) 



The proof of Theorem |3.4| is a simphfied version of the one given for Theorem |3.3| 



Only estimate (3.111 is needed in order to prove (i). 

4. Stationary solutions of the non local problem in a bounded interval. 

In this section we show the existence and uniqueness of a stationary solution of the 



non local problem (1.16) and analyze its relation with the travelling wave solution 
found in the previous section. We have divided the section in three subsections. In 
the first one we analyze the local problem (see (4.1)) in a bounded domain. The 



results from this subsection are used to obtain existence and uniqueness of stationary 
solutions for the non local problem (1.16). Finally, we show that these stationary 



solutions approach the travelling wave solution. 

4.1. The local problem in a bounded interval. In this subsection we study 



the existence, uniqueness and properties of stationary solutions of problem (2.4) when 



the domain is truncated to a finite interval. We impose non homogeneous Dirichlet 
boundary conditions, i.e., we consider the evolution problem 



Ut 



cu,. 



X e (a, b), t>0, 
u{a,t) = 0; u{b,t) = 1, t > 0, 
u(x, 0) — uq{x), X <E [a, 6], 



(4.1) 



with < Wq < 1 and analyze its set of equilibria lying between and 1. Observe that 
these equilibria are solutions of the boundary value problem. 



C/"(O + cC/'(O + /(C/(O)=0, a<C<6, 
U{a) = 0; U{b) = 1, < [/ < 1, 



(4.2) 



which can be interpreted as the first coordinate of the solution ([/, V) of the 2x2 
ODE 



U' = V, 

V ^~cV-f{U). 



(4.3) 



satisfying U{a) = 0, U{b) = 1, < U{£) < 1 for all ^ G {a,b), where we denote by 
- 

We will consider now several important properties of the solutions of the ODE 



(4.3). In the first place, we notice that two different solutions of (4.3) do not intersect 



due to the uniqueness of solutions of any initial value problem for this ODE. In what 
follows we will make use of this property without further notice. Observe also that 
since this ODE is independent of the "time" variable we have that if (C/(^), V{£,)) 
is a solution, then {U{^ + a), + a)) is also a solution for any a € JR. Hence, we 
may consider solutions of (4.3) defined for ^ € [0,r]. We will concentrate mainly 
on solutions starting at a point of the form {U,V) = (0,6*) and will analyze the 
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dependence and properties of this solutions with respect to 0, the constant c and 
so on. Since we are assuming that the nonhnearity / satisfies conditions (1.3), then 
system (4.3) has three distinguished solutions which are given by the stationary states 
(0,0), (1,0) and (a,0). 



There is another special kind of solutions of (4.3) which will play an important 
role in the analysis below and that we will denote as "simple solution" . 

Definition 4.1. A "simple solution" is a solution {Ui^),V{^)), ^ £ [0,r] of 
^n^ , joining (0,9) and (1,T) with 6*, T > (that is, ([/(0),y(0)) = (0,61) and 
(C/(r), y(r)) = (0, T) ) and with < t/(0 < 1, C e [0, r] (see Figure\t^. 

We start proving some results about these "simple solutions" . 

Lemma 4.2. // (U,V) is a "simple solution" then there exists vq > such that 

V > vq. Moreover, the curve {U{^),V{^)), ^ G [0,r], can be represented as a function 

V = P{U) forO<U <1 and P(C/) > for all < U < 1. 

Proof Observe that by definition it is V"(0) = >0 and V{r) = T > 0. 

If there exists a point ^~ e (0, r) with V{(,^) < 0, then let be the first point 
such that V{^) > for < ^ < ^0 and V{^o) = 0. By definition < r and V'{^o) < 0. 
Again we cannot have V'{^o) = since if this is the case, then V'{(,a) — — cy(^o) — 
f{U{£_o)) = —f{U{S,o)) — and this imphes that {U{^o),V{^o)) is an equilibrium 
point of ( |4.3| ) which is not the case. Therefore V'{£,o) < (hence f{U{£,o)) > 0) and 
in particular for ^ > but ^ near we have J7(^) < C/(^o) and V{^) < 0. Moreover, 
for < C < ^0 we have < U{^) < C/(^o) < 1 (if U{^o) = 1, since ^(^o) = we 
will have again another equilibrium point of (4.3)). Hence, the curve {U{£,),V{^)) for 
< C < Co joins the point (0,6^) with (C/(Co),0) with U{() increasing and V{^) > 0, 
that is this curve can be expressed as = P{U) for <U < t^(Co) for some function 
P. 

Consider now the region R ^ {(U^V) : < U < U(^o), -oo < V < PiU)}. 
Using the direction field of (4.3), we have that any initial condition starting in the 
set R either remains for all times in R or it leaves the region R through the boundary 
{(0, y) : — oo < < 0}. Since our solution satisfies U{b) = 1, necessarily has to 
leave the region R but then there will exist a point Ci e (Coj^) with C/(Ci) = and 
V{£^i) < 0, which implies that U{S,i + e) < 0, for e > small enough. This is a 
contradiction with the fact that < U{^) < 1. Hence V{^) > for all ^ G [0,r]. The 
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existence of the function P follows easily now from the fact that C/'(^) = V{£) > 
for aU ^ e [0,r]. □ 

Lemma 4.3. The following properties hold: 

(i) If V = P{U) is a simple solution then the time it takes to go from the point 
iUi^o),V{^o)) to {Ui^i),Vi^i)), that IS a - Co, is given by 



a -Co 



In particular the time it takes to go from {0,9) to (1,T) is given by 

r= 

(ii) If we have two "simple solutions" V — Po{U) and V — Pi{U) satisfying 
Pq{U*) < Pi{U*) for some U* G [0, 1], then Po(U) < Pi{U) for allO<U <1. 

(Hi) If we consider ( |4.3[ ) for two constants cq < c\ and simple solutions V = Pi{U) 
for the system with Ci, i = 0,1, respectively, then if Po{U*) = Pi{U*) for some 
< U* < 1, it is Pq{U*) > P{{U*). In particular the orbits can only cross once (see 



Figure 4-2) 



(iv) If > \c\ + \ f\oo + 1 where \ f\oo — max{|/(s)|, < s < 1} then the solution 
starting at (0, 9) is a simple solution joining (0, 9) with (1, T) for some T > — |c| — 
|/|oo- Moreover this simple solution V — P{U) satisfies P{U) > 9 — (|c| + |/|oo)C^ for 
< U < 1 . In particular 

U, as 9 +00. 



PiU) 

(v) There exists 9fj>Q such that the solution starting at (0, 9) for each 9 > 9q is 
a simple solution V = Pe{U) with 

du I 
-^^^^+00, as 9 ^9^. 

(vi) For any 9q > there exists a cq = co(0o) > such that for all c > cq any 
"simple solution" starting at {0,9) with 9 > 9o ends at (1,T) with T < 9. Similarly, 
there exists a ci = ci{9o) < such that for all c < ci any "simple solution" starting 
at (0, 9) with 9 > 9o ends at (1, T) with T > 9 



Proof, {i) It follows from standard results. Notice that the first equation of (4.3) 



reads ^ = V ^ P{U) which implies = 'pfu) ^o^^ where we obtain the result via 
integration. 

{ii) This follows from the uniqueness of solutions of the ODE ( |4.3[ ). 

{Hi) Note that = -c + f{U)/P{U). Hence, if Po(C^*) = Pi{U*), then 

P(,{U*) = -co + f{U*)/P{U*) > -ci + f{U*)/P{U*) = P[{U*) 

{iv) In the {U,V) plane, the segment joining {0,9) with (1,T*) is given by = 
{r*-9)U + 9, for < [/ < 1. Hence, if 6* > |c| + |/|oo + 1 and T* = 61 - |c| - |/|oo, the 
straight line is y = — (|c| + \ f\oo)U + 9 and notice that we have F > 1 for all {U, V) 
in the segment. But, in this segment, the direction field of the ODE ( |4.3| ) points 
"toward the right" of the straight line, since the scalar product of its normal vector, 
(|c| + I/loo, 1), with the direction field {V, -cV + f{U)) is 

(|c| + \f\oo)V -cV + f{U) - V{\c\ -c) + \f\ooV + f{U) > 
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V 




(0,0) 



V=Pg(U) 



u 



V=P^(U} 



(1,0) 



Fig. 4.2. Two simple solutions V = Po{U), V = Pi{U) for cq < ci which crossing. 



(|c| + |f| , 1) 
an 



v=-(|c|+|f| )u+e 




(0,0) 



(1.0) 



Fig. 4.3. The segment and the direction field. 



since if (U, V) lies in the straightline then V > I, see Figure 4.3 



This imphes that if we start at (0, 9) then the orbit of the ODE cannot cross this 
hne segment and therefore we wiU have P{U) > 9 — (|c| + |/|oo)C^ for a\\ < U < 1 
and in particular it will arrive at a point (1, T) with T > T*. The last statement of 
this item is obvious from the bounds of P{U) and taking 9 — > +oo. 

(v) It is clear from {iv) that there exists > such that the orbit starting at (0, 9) 
for 9 > 9i are all simple solutions. Moreover, from Lemma |4.2| and the continuous 
dependence of the solutions of the ODE with respect to the initial conditions, we have 
that if the solution starting at a point (0,9*) is a simple solution, then there exists a 
small e > such that for each 9 E {9* — e, 9* + e) the solution is also a simple solution. 
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(1,0) 



Fig. 4.4. Phase planes for Lemma\4--4 



Let US define 



Of) = 'mi{6 : the orbit starting at (0, O^yiO^ > is a simple solution } 

Then, by definition the orbit starting at (0, ^o) is not a simple solution. But, neces- 
sarily this orbit satisfies < C/ < 1 and, if it reaches a point of the form (1, T) in finite 



time, by Lemma 4.2 it will also be a simple solution. Therefore by the continuity of 
the fiow, we necessarily have either or T = 0, which leads to 



du 



+00, as 9 ^ 9i 



(vi) Observe that if > and > 9q then, on the horizontal line V — 9, the 
vector field of the ODE is (9, -c9 + f{U)). Therefore, if we choose cq = l/U/^'o > 0, 
then for any c> cq the vector field "points downwards" in this line. This implies that 
if we start with a solution at (0, 9) and the solution travels to a point (1, T) then the 
whole orbit lies below the straight line V — 9. Hence, T < 9. Similarly for the other 
case. □ 

We also have another two distinguised solutions with < ?7 < 1 and V > 0, 
which are usually denoted as the "unstable orbit" from the equilibrium (0, 0) and the 
"stable orbit" to the equilibrium (1,0). These are special solutions since they are 
defined either in an interval of the form ^ G (~oo, rg) or ^ S (ri, +00). 

Lemma 4.4. There exists a unique value c* € M such that the following holds 

(i) If c < c* the unique "unstable orbit" emanating from (0,0) with < U < 1 
ends at a point (Uo,0) for some value < C/q < 1 and the unique "stable orbit" 
converging to (1, 0) with < U < 1 starts at a point of the form (0, ^o); with 9q > 0. 

(ii) If c > c* the unique "unstable orbit" emanating from (0,0) with < U < 1 
ends at a point (1, Ti) for some value Ti > and the unique "stable orbit" converging 
to (1, 0) with < U < 1 starts at a point of the form [Uq, 0). 

(Hi) If c = c* , both orbits coincide and we have a unique orbit with < U < I and 
V > which comes out from (0, 0) (as —00 ) and comes in to (1, 0) (as ^ +00 ). 

Proof. The proof is a simple exercise of phase plane techniques and we leave if 
for the reader. See Figure 4.4 and the book [H] for details. □ 

Remark 4. Observe that the value c* from the previous Lemma is the speed of 



propragation of the travelling wave of problem (1.1 1. 
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With the resuhs above we can prove the first result on existence of solution of 



(4.2) with c fixed. 

Proposition 4.5. For every c G M, there exists a unique solution <i>c of (4.2 1. 
Furthermore, $c is strictly monotone, that is, ^'(x) > for all x G [a, &]. 

Proof. Observe that $c is a solution of ( |4.2[ ) if and only if ([/, V) = (<i>c, $c) is a 
simple solution as defined above. But if we consider the value from Lemma 4.3 {v) 



and the simple solution V — Pe{U) which starts at (0, 9) then the time it takes to go 
from (0,F(0)) to (1,P(1)) is given by 



r{B) = 



dU 

mj~y 



which by Lemma 4.3 [iv) and (w) satisfy that there exists at least one 9 for which 
r{6) = b — a. From Lemma 4.3 (ii) we have the uniqueness. □ 

As we show below, the monotonicity of the equilibrium, yields its asymptotic 
stability. 

Proposition 4.6. Let c e M and the equilibrium solution $c of (4.1). Then, 
there exists K,a > such that for \\uq — $c||oo small enough it holds 



\ui-,t) - $c||oo < Ke' 



(4.4) 



for u the solution of (4.1). 



Proof. We consider the linearization about $c of (4.10) 



Ut = + CUx + f'{<^c)u, 

u{a, t) — u{b, t) = 0. 



We need information about the spectrum of the operator 

Lo<Z:=g"+c<z' + /'(<fc)(?, 



(4.5) 



(4.6) 



in D{Lo) = {(? e C^[a,b] : q{a) = q{b) = 0}. Deriving the equation in (4.2 ) we obtain 
that (/) = <i>J, > satisfies the boundary value problem 



(j)'{a)+c(j>{a)^{), (l)'{b) + c(j){b)=0, 



(4.7) 



where the boundary conditions are obtained by evaluating ( |4.2[ ) at ^ = a and ~ b 
and using that /(O) = /(I) = 0. The change of variables z(^) — e^^q{£,) in (4.7) leads 
to the self-adjoint problem in L^(a, b) 



-/'($,) z = 0, 



(4.8) 



z'(a) + -z(a) = 0, z'{b) 



-z{h) = 0, 



For (4.8), the positive mapping ^ — >■ ei^(t){C) > is an eigenfunction associated to the 
eigenvalue 0. Then, by Krein-Rutman's Theorem, is the smallest eigenvalue of the 
operator 



Bq := -z" 



-/'(^JU, 



(4.9) 



20 



for the boundary conditions in (4.8). Then, 



mm 

_f/i(a,b) 



< mm 

Hiia,b) 



SO that the smaUest eigenvalue of B with homogeneous Dirichlet boundary condition 
is strictly positive. This implies (T(io) C {z G C ; Rez < 0}, yielding the asymptotic 
stability of $c. □ 
Corollary 4.7. With the notations above, we have (7{Lq) C {z G C : Re z < 0}. 

4.2. The nonlocal problem in a bounded interval. In this section we study 
the existence of equilibria of the non local equations (1.9) and ( |1.15 ) restricted to a 
finite interval with nonhomogencus Dirichlet boundary conditions. It is immediate to 
see that the resulting Initial and Boundary Value Problems (IBVP) for both equations 
are the same, namely: 



Mz{-)) 



+ X{ux)ux + f{u), 
F{1) 



X e (a, 6), t > 0, 



l^(-)ll^ 



L2(a,b) 

u{a,t) = 0; u{b,t) = 1, t > 0, 
^ u{x,0) ~ uq{x), X e [a,b], 



(4.10) 



We have the following 

Proposition 4.8. Problem (4.10) is locally well posed, in the sense that for any 
initial data uq € H^(a,b) with uo(a) = 0, wo(fe) = 1 there exists aT = T{uo) and a 
unique classical solution u{x,t) defined for the time interval [0,T(uo)). 

Moreover, if < uq < 1, then the solution u(x,t) is globally defined and it also 
satisfies < u{x,t) < 1. 



Proof. Observe that problem (4.10) can be rewritten as a more standard problem 



with homogeneous Dirichlet boundary conditions with the change of variables: u{x) = 
w{x) + h{x), where the function h{x) = |5f ■ Notice that if u{a) — 0, u{b) — 1, then 
w{a) = w{b) = and problem (4.10) takes the form 



Wt^Wxx + F{vu) X € {a,b), t > 0, 

w{a,t) = w{b,t) =0, t > 0, 

w{x, 0) = wq{x) = uo{x) + h(x), X e [a, b]. 



(4.11) 



where F : Hl[a, b) — > LF'[a, b) is the map defined by F[w) — \{vux + j;z:^){wx + ^^) + 
f{w + h{-)). We can easily see that this map is well defined, since if w G -ffg(a, 6) 
evaluating in the following expression, we have 



1 



b-. 



||2 11 



L2(a,6) 

and therefore the denominator in the function A is bounded away from and \{wx + 
j3^) is well defined. Moreover, following standard arguments, the map F is Lipschitz 
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Fig. 4.5. The equilibrium in a bounded domain: (Left) In the phase space {U, V); (Right) As a 
function of x . 



on bounded sets of Hl{a, b) which from standard techniques, see [TOl, we obtain that 



problem (4.11 ) is locaUy well posed in i?o(a, h) and therefore problem (4.10) is locally 
well posed for any initial condition ug G H^(^a,b) satisfying uo{a) = 0, uo{b) — 1. 
Standard regularity results applied to ( 4.10[ ) (notice that A is independent of x) show 



that the solution is classical for t > 

If we consider now that < uq < 1, then we may argue by comparison with the 
constants to show that as long as the solution exists, it will also satisfy < w < 1. 
Let us argue by contradiction. Assume the solution is negative at some time T > 0. 
Then, for e > small enough there exists a < < T such that u{x, t) > —e for all 
X G [a,b], < t < and there exists x^ G (a, b) such that u{x^,t^) — — e. This imphes 
that ut{xe,t^) < 0, Uxxixe^te) > 0, Ux{xe,te) = and f{u{xe,t^)) = /(-e) > 0. 
Which is a contradiction. In a similar way we may proceed with the upper bound 
u < 1. This shows that, as long as the solution exists we have < u{x,t) < 1. 
With standard continuity arguments we show that the solution is globally defined 
and satisfy the bounds. This concludes the proof of the proposition. □ 

Remark 5. Notice that although we have been able to obtain a comparison result 
with the constants and 1, we do not have a general comparison argument for equation 



(4.10). That is, if the initial conditions are ordered uq < vq we cannot conclude that 
the solutions satisfy the same ordering for positive times. The lack of these comparison 
arguments for this equation is very much related to the lack of maximum principles 
for the associated linear non local operators. This lack is a serious drawback, specially 
when analyzing the stability properties of the equilibrium of this equation, see Remark 
below. 

The following result provides a characterization of the stationary solutions of 



(4.10) 



Lemma 4.9. A function $(•) e H^{a, b) with $(a) = 0, $(&) = 1 and < <I>(x) < 



1 is an stationary solution of (4.10) if and only if $(•) is a solution of (4.2) with 
c — A($') and satisfies 

$'(a) = $'(6). (4.12) 



Proof. It is clear that any equilibrium solution $ of (4.10) satisfies (4.2) for 
: A($2;). Multiplication by $' in the ODE in (4.2) and integration along (a, 6) 
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yields 



Using that (/($),$') = F{1) and the definition of A in (4.10) we obtain that $ 
satisfies ( |4.12| ). □ 



We can proceed now to prove an existence and uniqueness result for stationary 



solutions of (4.101 



Theorem 4.10. There exists one and only one stationary solution $ of (4.10) 
with < < 1. Moreover, this solution is strictly monotone increasing in x. 



Proof. From Lemma 4.9 a stationary solution $ of (4.10) with < < 1 is 

(4.13) 



the first coordinate of a "simple solution" of the ODE 



V' = -cV-f{U) 

where c = A($') = -^t^^^ — G M. Lemma|4.2|proves that $'(a;) > for all x e [a, b] 



and therefore $ is strictly monotone increasing. 

Uniqueness is obtained as follows. Assume that there exist two solutions $i, $2 
and denote by ci — A($i), C2 = A($2)- From the uniqueness of solutions given in 



Proposition 4.5 we get that ci 7^ C2. Then, if we denote by {Ui,P{Ui)), {U2, P{U2)) 
the two simple solutions associated to $1 and $2 respectively, from Lemma [4. 3| (Hi) 
and from Pi(0) = Pi(l), ^2(0) = P2(l) we must have that either Pi{U) > P2{U) or 
Pi{U) < P2{U) for all < ?7 < 1. In both cases we have 

dU ^ dU 



Pi{U) Jo P2{U) 

which is a contradiction. This shows uniqueness. 

Existence is shown as follows. We know from Proposition |4.5| that for every fixed 
Co and r = b — a there exists a unique solution $co of ( |4.2[ ) , which actually is given by 
$co = Uq where {Uq,Vq) is a "simple solution" joining (0, 0o) with (l,To) for some 
Bq^Tq > 0. If it happens that 6*0 = Tq, that is ^'cgio-) = ^'cgi^)^ then Lemma 



4.9 



shows that this function $^0 is the stationary solution we are looking for. If 6*0 7^ 1 0, 
let us assume that 9q < Tq (the other case is treated similarly). For c > cq and the 
same r = b — a, again by Proposition |4.5| we have the existence of a solution $c which 
again is given hy = U where (U, V) is a "simple solution" joining (0, 9) with (1, T). 
But since r is the same for both solutions and we have r — pfiu) ~ fo p (U) then 



necessarily, both solutions must cross at least at some point and by Lemma 4.3 they 
can only cross at one point and it must be satisfied 9 > 9^^ T < Tq. Moreover, from 
Lemma [4?3| we can choose Ci > cq large enough such that for this value ci the unique 
simple solution joining a point of the form (0, 9i) with (1, Ti) in a time r = b — a 
satisfies > Ti. By the continuous dependence of the solutions $c with respect to 
the parameter c, we will have that there will exist a value c* € (co,ci) such that the 
unique solution $c* travelling for a time r = b — a joins a point of the form (0,9*) 
with (0, T*) with 9* = T* > 0, that is $'(0) = $'(fe). This is the desired solution. □ 

4.3. Convergence of the stationary solutions to the travelhng wave as 
the length of the interval goes to +00. In this section we will pass to the limit as 
the interval grows to cover the whole line and we analyze how the solution encountered 
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in Theorem |4.10| behaves as the length of the interval goes to infinity. The first step 
is to prove the convergence of the wave speed to the one of the travelling wave. More 
precisely, 

Lemma 4.11. Let Aj. be the unique value given by Theorem \4.10\ for which an 



equilibrium of (4.10l-(4.12| exists on the interval {a, b) with r = b ~ a. Then, 

- Aool 0, as r^+oo. (4-14) 
where Aoo — c* from Lemma \4-4\ ^hat is, the speed of propagation of the travelling 



wave of equation (1.1 1 



Proof. Observe first that the value of A^ really depends only on r = 6 — a and not 
on a or b. 

Assume that the result is not true. Then, there is a sequence {r„}„gp} with 
r„ — )> cxD, as n — >■ cx), and e > so that if we denote by A„ :— Ar„, then either 
> Aoo + e or A„ < Aqo — e. So let us assume that A„ > Aqo + £, for all n, the other 
case is treated similarly. 

Observe that from Lemma |4.4[ in the phase plane associated to the equation 



(4.10) for c = Aoo + e there is an orbit ($*,$'^) arriving at (1,0) from (0,T*), for a 



certain T* > 0. This orbit is also represented as V — P*(U) for < U < I. By (iii) 



of Lemma 4.3 and (4.12), none of the orbits ($„, $^), which are given hy V — Pn{U), 
can cross V — P*{U) and it has to be Pri(O) = T„ > T*, for all n. It follows then that 
Pn{U) > P*{U) for aU < J7 < 1. Furthermore, the graph of the function V = P„(C/) 
is also above the straight line passing through (1, T*) with slope (1, —Aoo + e)- This 



comes from the fact that for a in (1.3 1 and a < u < 1, it is f{u) > and the field 
in the phase plane is proportional to (1,— A — f{u)/v) with —A — f{u)/v < —A. It 
follows that the orbit V = Pn{U) has to arrive at (1,T„) from above this line. But 
then, the time it takes to travel from (0, T„) to (1, T„) remains bounded, i.e., by (i) 
of Lemma 14.31 it holds 

Jo ^n{U) 

This is in contradiction with the fact that r„ — > +oo. □ 
Lemma 4.12. Let $r be the equilibrium obtained in Theorem \4.lO\ in the interval 

(0,r). Then the orbit (^ri^r) *^ the phase plane converges to the orbit associated to 

the travelling wave on the whole line, (<I>oo,$^) as r —> oo. 

Proof. Observe that the orbit (<I>r,$^) is a simple solution and it is given as 

V — Pr{U), < U < 1. Moreover, we know that the travelling wave ($oo,*5j^) is 

given as the function V = Poo{U) for < t7 < 1. We will show that Pr — Poo as 

r —7- +00. 

Assume the lemma is not true. Then we will have a sequence of r„ — +00 and a 
Uo e [0, 1] such that Fr„ {Uq) Vq > Poc,iUo)+S for some S > 0. Notice that we have 



used the fact that Pr > Poo- But we know from Lemma 4.11 that Ar„ — J> Aoo- Hence, 
by continuous dependence with respect to the initial conditions and with respect to 
the parameters appearing in the equation, the orbit ($r„ , ) converges to the orbit 
of the ODE with A = Aoo passing by {Uq, Vq). Since Vq > Poo{Uo) + S we have that 
this orbit takes a finite time to go from the line J7 = to the line U — 1. This is a 
contradiction with the fact that r„ — > +00. □ 

We will normalize the orbit ^r) so that the time ^ = will correspond to the 
unique point for which $r(0) = 1/2. Hence, we will denote by a{r) < < b{r) so that 
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Fig. 4.6. Convergence of the equilibrium to the travelling uiave 



b{r) - a{r) = r and ($,.(a(r)), $;(a(r))) = (0, 9) and KiKr))) = (1, 0). In 

a similar way we may normalize the travelling wave solution so thar $oo(0) = 1/2. 
We have the following 

Proposition 4.13. With the notations above, we have both, 
a{r) — > — oo, and b{r) — > +oo. 



Proof. Assume one of them is not true. For instance, let us consider that there 
exists a sequence r„ — >■ +oo such that 6(r„) bg < oo. This implies that the finite 
interval [0, b{rn)) approaches the finite interval [0, to) and therefore by the continuous 
dependence of the solutions of the ODE with respect to the parameters and the initial 
conditions in a finite time interval [T^, we will have that ($r„ ^'rn (^(''n))) = 

(1,T„) — > (<I>oo(&o), ^^x)(^o)) and this implies that <i>oo(&o) = 1, which is impossible 
for any b^ < oo, since $oo is the travelling wave solution. 

A similar proof shows that a(r) — > — oo. □ 

We may also prove 

Lemma 4.14. With the notations above, if we extend the function {$.r} by to 
the left of a{r) and by 1 to the right of b[r) (and we still denote this function by 
then 

||$r - $oo|ki.~(E) + \\K - ^LWl^Vl) ^0, as r ^ OO. 



Proof. The convergence in follows directly from Lemma 

over, notice that since A,- — >■ Aoo and using that A,. = — F(l)/||$^||^2(j 



4.12 



More- 



and A^ 



-J^(l)/ll*'oolli2(R3 we have that 



i*:>oiii.( 



Hence, consider a small enough parameter e > and let us fix a large enough 



interval [-T,T] such that ||$^|||2( 



-T,T)) 
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< e. Then, from the convergence of the 



orbits given by Lemma 4.12 we have that hm^-j-oo ll'i'r ^ ^'oo\\l'^{-tt) ~ which 



imphes that Umr_j.oo \\^'r\\L^{-T,T) = II*^^IIl2(_7-,t)- Hence, 

Urn \\K\\h(s.\(~T.T)) = lim \\K\\h(R)~ 1™ \\K\\h(~T.T) 

and therefore, 

hm ii<i>;-$:,oiii2(R) < lim iii>;-<i>:,oiii2(_T,T) 

+2,1™ \\K\\h(R\i-T,T)) + '^WooWIhrm^t^t)) < 4e- 
Since e is arbitrarily smah, we show the Lemma. □ 

5. Asymptotic stability of the stationary solutions of the nonlocal prob- 
lem. We analyze in this section the stability properties of the unique stationary 
solution of the nonlocal problem (4.101 in the bounded domain (a(r),6(r)). We con- 
sider the normalization of this equilibrium explained in the previous subsection, that 
is <l>r(0) = 1/2 and to simplify the notation we will denote the interval by (a, 6) in- 
stead of (a(r), b{r)), unless it is necessary to specify the dependence of the domain in 
r. 

The linearization of ( |4.10 1 around is given by, 

Wt=W^^ + \{<i>r)w^ + f'{^r)w+Trr{w)^'. xe{a,b), t>0, , , 

wia,t) = w{b,t) = 0, ^ > 



with TT,. the linear nonlocal operator 



^.H = -2A($.)^^j|^. (5.2) 



The equilibrium $r will be asymptotically stable if the spectrum of the linear 
operator L'' : D{U) C L'^{a,b) L'^{a,b) wilYi D{U) = H^{a,b)n H^{a,b), given by 

L''w; := Wx.j; + X{^r)uix + f{^r)w + nr{w)^'r (5.3) 



is contained in the left half of the complex plane. We recall that, by Proposition |4.6[ 
this is the case for the local operator 

L^w ::^w^^ + X{<Pr)w^ + f'{^r)w, w€H^{a,b). (5.4) 

Observe that L'^w ~ L'QW+Trr{w)^'^ and the operator w — Trr{w)^[. has 1-dimensional 
rank and can be expressed as 

2A($^ f 

This operator is of the form w — > A{w,B) with A{-) = and B{-) — 

II '^'r II 2 

<!>"(•) = -A(<i>r)$r(') ~ /(*i-(-)) a^^d it is a bounded operator from to with 
finite rank. Several properties of the operator are inherited from the operator Lq-. 
both operators have the same domain, both operators have compact resolvent and 
therefore the spectrum is only discrete, formed by eigenvalues with finite multiplicity. 
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Nevertheless, all the eigenvalues of operator Lq are real (there is a standard change of 
variables transforming Lq to a selfadjoint operator) but the operator L*" may not have 
this property. Actually, unless A = B operator is not selfadjoint. There are several 
studies of the spectrum of operators of the form w — !■ Lq{w) + A{w, B) but none of 
them guarantee us that for our particular case, the spectrum lies in the half complex 
plane with negative real part. Actually, with the known results in the literature we are 
not even able to show that the spectrum of L is real. See [HI HOI [HI [3 [E] for results 
in this direction. One important observation is that in the case that the interval is 
the complete real line, that is r = oo, then is the eigenfunction associated to the 
eigenvalue for the operator and therefore the operator = L"^ + 7roo(w)$^ 
has an special structure that will allow us to show that (t(L°°) — a{L'^) and that 
G a{L°°) with multiplicity 1. As a matter of fact this will give us an alternative proof 
of the asymptotic stability (with asymptotic phase) of the travelling wave solution of 



the nonlocal equation in the whole real line (see Theorem 3.3 1. The fact that <i>J, is 
not an eigenfunction of Lg, for finite r (as a matter of fact <1>J, does not even satisfy 
homogeneous Dirichlet boundary conditions) will not permit us to perform a similar 
argument in a bounded interval. Paradoxically, the analysis in the whole real line is 
"simpler" than the analysis in a bounded interval. 

Nevertheless we will be able to prove the asymptotic stability of the stationary 



solution of the non local problem (4.10) for large enough intervals using a perturba- 



tive method. The proof is divided into three parts. In the first one we prove some 



properties of the spectrum of the non local operator (5.3) on the finite interval (a, 6). 
We next fully analyze the spectrum of the limit operator on the whole line M. Finally, 
we prove the convergence of the spectrum of L*" to the spectrum of L°° as r — > +oo. 

5.1. Spectral properties for any finite interval. The results in this section 



apply to the stationary solution obtained in Theorem 4.10 in the finite interval 
(a,&). 

Let us start with a general and rough estimate of the spectrum of L*" but which 
is uniform for all r > 1. 

Proposition 5.1. There exist po G M+ and cj) e (tt, 27r) such that if we define 
the sector Sp(,,0 = {z £ C, \Arg{z — po)\ > cp}, then (j{L^) C Sp^^^ for all r > 1. 

Proof. Note that ^ £ <^{L^) if and only if there exists u G H^{a,b) n Hl{a,h) 
such that U'u — fiu. But the operator L'^ can be written as L'^u = Au + N{u) where 
N : iJo(a, &) L'^{a,b) is defined as Nr{u) = \{^r)ux + f'i^r)u + 7rr(u)$J, and as 
usual Am — Uxx- Observe that from Lemmas [4 . 1 1 1 and [4 . 1 4| we have that the operator 
Nr is bounded uniformly in r for r > 1, that is, there exists a constant Cq independent 
of r = b— a such that || A^ru||L2(Q < Co||u||//i(a,6)7 for all r > 1. 

On the other hand, standard estimates using the spectral decomposition of —A 
with Dirichlet boundary conditions in (a, b) show that for fi ^ 

\\i-^ + liirTc(LHaM,HU-M) ^ dist(^,M-) + + 
Hence, fixing (f> G {tt, 2tt) we can choose po > large enough so that we have 
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|(-A + /i/)- 



\\c{L^ia,b),Hai'>;b)) - (2Co)2 ' ^ ^''O''^' 



Therefore, if /i G C\SpQ^0 and if there exists u £ H^{a, b)r\H^{a, b) such that L^{u) = 
/iu, then, u = N o (—A + pT)^^u which implies that ||ti||i2 < ||-/V||£(//i_l2) ||(— A + 
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/^-^) "'^ll£(L2.ifi)||'t^||L2 < C'o2g^ll'"lli^ — and therefore u = 0, which imphes 

that n ^ (t(l°). ° □ 

This rough estimate of the spectrum of L"^ ahows us to prove that if there is 
an eigenvalue of L"^ with positive real part, then necessarily we will have that it is 
uniformly bounded in r, that is 

Corollary 5.2. With the notations of the previous proposition, we have that 
for any value > 0, we have 

{z e cr{L''), Rez > -J/} C {z e C, -i^ < Rez < po, \Im{z)\ < (po + i^) sin(0)}. 



Lemma 5.3. Let fi e a(L^) p{Lq) . Then, fi is at most a geometrically simple 
eigenvalue of , that is, Ker{L^ — pi) is one dimensional. Moreover, the associated 
eigenspace is generated by y, the unique solution of 

(LS-A^/)y-<, (5.5) 



(v' 

^r{y) = -2X^j^^-l. (5.6) 
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Proof. Let w 7^ be such that L''w — pw. Then, 

- pl)w = {LI - pl)w + 7r^(w)$; = {LI ~ pl){w + 7r^(u;)y), 

so that 

w = —Trr{w)y. (5.7) 

The above implies that p is at most a simple eigenvalue of L^ with eigenspace gen- 
erated by y. Identity (5.6) follows after applying the linear operator tt^ in (5.7). 
□ 

However, it will be still useful for the last part of our argument. 

Proposition 5.4. There is no real eigenvalue p > in a{L^). 

Proof. Let us assume that there exists an eigenvalue ^ > of L''. Since cr(Lo) C 
{z g C : Rez < 0} (see Proposition 4.6 and Corollary 4.7) then p e p{Lq). Let y be 
as in (5.5) for this value of p. Then, by definition, y satisfies 

y" + Xy' + {r{^r)-t^)y^K- (5-8) 
Multiplying in ( |5.8[ ) by <i>J, and integrating in (a, b) we have 

(y", $') + X{y', $;) + (/'(<i>.)y, K) - i^iy, K) = Wrf- 



But 



so that 



{f{<^r)y, K) = ~{y', /(<f r)) = (y', K + AO 



(y",0 + (y',<i>;')+2A(y',$;)-^(y,ci>;) 
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1$: 



' l|2 



Then, by (5.6), it holds 



(y",<i>;) + (y',0 =M(y,$;> <0, (5.9) 
where the inequahty foUows from the maximum principle applied to --Lgy+i-iy — 



and taking into account that <I>J, > 0, so that y < in (a, 6). But, from Lemma 4.9 



we know that 9 $J,(a) = > 0, which implies, together with (5.9), that 

(y", <) + (y', O = y'Kt = 0{y'{b) - y'{a)) < 0, 

and therefore y'{b) < y'{a). But, on the other hand, the fact that y < in (a, 6) 
together with y = in a; = a, 6, imply that y'(a) < < y'(b). Therefore y'{a) = 
y'{b) — 0. But this is impossible, since if, for instance, y'{a) = 0, then y is a solution 
of the initial value problem LqJ/ ~ "I>J, in (a, b) with y{a) — y' {a) = and this implies 
that y"{a) = ^r(a) > 0, so that with x near a we have y > which is not true. □ 

Remark 6. i) This proposition would be enough to finish the proof of the asymp- 
totic stability if the non local operator U" had the property that the eigenvalue with 
the largest real part were real. For instance, this could be obtained if U" satisfies the 
hypothesis for a Krein-Rutmann type of theorem. But for this theorem we need to have 
maximum principles and are unable to prove this principles for this nonlocal operator. 

a) Observe that this proposition does not exclude the possibility of having complex 
eigenvalues with positive real part. Actually, we will be able to exclude this possibility 
only for large enough intervals by using a peturbative argument. The fact that this 
operator may present complex eigenvalues with positive real parts for some intervals 
(a, b) is an open interesting question. 

5.2. Spectrum of the nonlocal problem in the whole line. In this sec- 
tion we analyze in detail the spectrum of the corresponding nonlocal operator in the 
entire real line. This operator is the one associated to the linearization around the 
asymptotic equilibrium $oo, that is, 

L^{w) = w^^ + A($oo)wx + f{^oo)w + ^ooH^^x), (5.10) 
where now t:^ stands for the linear operator 

-2A($ 

2 

We will use several important properties of the spectrum of the local operator 

L'^w.^w" + \{^^)w' + f'{^^)w, (5.12) 

We have the following. 

Lemma 5.5. With respect to the spectrum of , defined by (5.12), we have 
(i) The essential spectrum aess{L'^) C {z G C : Rez < max{/'(0), /'(I)}}. 
(a) There exists Q < v < - max{/'(0), /'(I)}} such that a{L^) n {z € C : Rez > 
—v} ~ {0} and the eigenfunction associated to fj, = is 

(Hi) There is no solution w S -ff^(M) of L'^w — ^'rx>- Therefore, is an algebraically 
simple eigenvalue of Lf^ , that is, Ker[{L^Y) = Ker{L'^) =span{^'^} 
We refer to Appendix [B] for a proof of this result. 

Both and L°° are sectorial operators and are related by 

L°°{u) = L^{u) + n^{u)-<^'^ (5.13) 
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^oo{w) = —^j^{U\w^). (5.11) 



In the following Proposition we show that L°° enjoys the same spectral properties 



as 



Proposition 5.6. Let L°° be the linear operator defined above in (5.131. Then 
a{L°°) = ct(L3°) and G a{L°°) is an algebraically simple eigenvalue of L°° . In 
particular the three items (i), (ii) and (Hi) from Lemma \5.^ also hold for L°° . 

Proof By applying integration by parts it is easy to see that tToo{^'oo) = 0- Then, 
since L^^'^ — 0, we also have that i°°<I>^ = 0, so that G cr{L°°) with associated 
eigenfunction . In order to see that is a simple eigenvalue of , let us consider 
first (f> G D{L°°) with L°°(j) = 0. In case 'noo{(f>) = 0, then it is = and, since is 
a simple eigenvalue for ig°, it must be </> ~ Let us assume now that 7roo((/>) 7^ 0. 
Then, it follows ijf </) — ~^Too{<f)^'oo^ which is impossible by Lemma 5.5 [Hi). With 
a very similar argument it is possible to show that there is no w G satisfying 
L'^w = (j)'^. Hence, is an algebraically simple eigenvalue of L°° . 

We show now that p{L^) C p{L°°). So, let ^ e p{L^), / G X and wj e D{L^) 
be the unique element of D{L^) such that 

L^Wf - pwf = /. (5.14) 

If T^ooiwf) = 0, then L°°Wf — pwf — f . In case TTooiwf) 7^ 0, we can consider 

w* ^Wf + -n^{ws)<^'^, (5.15) 

since we already know that /i 7^ 0. Then, using that ig°<I>^ = 7roo(^'^) — and 
(5.13), one gets 

L°°w* - pw* = f. (5.16) 

This proves that L°° — jil is onto. 

Let us assume now that there exist two elements w\,W2 G D{L°°) with 

L^w*-iiw*=f, for j = l,2. 

Then 

L'^w* ~ pw* = -7roo(u)*)<I>^ + /, for j = 1, 2. 

From the above it is clear that tTqo {w^ ) = tTqo (wj ) implies — since — fil 
is one to one. In case TTooiwl) ^ iToo{w2), we consider w — wl — W2 and we get 
L^w - fiw = -7roo(w)$J^,, which implies w = -7roo(w)(Lg° - fj,I)~^^'^. But (Lg° - 
/Lt/)~^<i>J^ = —^'^/fj, and therefore w ~ and 7roo(w) = 0, which is a contradiction. 
The fact that {L°° — p,I)^^ is bounded is clear from the expression 

(L- - /i/)-l/ = (L- - + -^eo(ii'/)$:x,- 



which is obtained from (5.15)-(5.16). This shows that p{L'^) C p{L°°). 

The proof of the other inclusion p{L°°) C p{L'^) is completely symmetrical to 
this one, once we know that is also a simple eigenvalue of L°° . We just need to 
express L[f = L — tToo'I'^ and remake the proof we have just shown. □ 

Remark 7. (i) In particular, we have that the spectrum of L°° apart from is 
located in the left half plane, i.e., there exists v > 0, such that 

sup{Rc^ : ^ G (T(i°°), p^O} = -iy, for some > 0. (5-17) 

(ii) Observe also that since the operator u ■Kao{u)^'r^ is a compact operator 
then aess{Loo) = fTes.(iSo) C {z G C : Rez < niax{/'(0), /'(s)} }, see 
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5.3. Spectral convergence and asymptotic stability of the stationary 
solution. In this subsection we will end up proving the asymptotic stability of the 
stationary solution of the non local problem. We will obtain this via a convergence 
of the spectrum of the operator L'^ to L°° . In order to prove this spectral convergence 
we will use the theory of regular convergence developed in [29l |30l [31] and the related 
results in [4]. The necessary definitions are given below. 

Let E and F denote separable Banach spaces and let {£'r}r>o and {Fr}r>a be 
families of separable Banach spaces. Let {pr}r>Oi Pr G ^{E,Er) and {(?r}reN, 9r G 
C{F, Fr) be linear bounded operators such that 

limr_>.oo Ibrellfi^ ||e||B, for every eeE and 
lim^^oo Ikr/llF, WIWf, for every / G F. 

A family {er}r>Oi G E,., is said to be 7^ -convergent to e G -E, written — ^ e, if 

lim \\er — Pr^WEr — ^- (5.19) 

A family {er}r>0: &r G £'r, is said to be -compact if every infinite sequence con- 
tains a -convergent subsequence. Analogous definitions apply for Q -convergence 
and Q -compactness. 

A family of bounded linear operators {Ar}r>Oi € £{Er,Fr), is said to be 
T'Q -convergent to A € C{E,F), written Ar — > A, as r ^ oo, if — >■ e implies 
Ae, as r — > oo. The -convergence is said to be regular if for every bounded 
sequence {er}r>Oi G Er, such that the sequence {Arer}r>o is Q-compact, it turns 
out that {er}r>o is 7^ -compact. 

The relevance of the VQ regular convergence is that we obtain the following 
result, which is taken from [29l|30l[3l] in a simplified version. 

Theorem 5.7. Assume we have the family of operators A{s) — A — sB G C{E, F) 
and Ar{s) — A^ — sB^. G C{Ej.,Fj.) where the parameter s £ S, a bounded subset of 
the complex plane C, which satisfy the following hypotheses: 
(i) Aj.{s) VQ -converges regularly to A{s) for all s € S. 
(a) For each s € S the operators Ar{s) and A{s) are Fredholm with index 0. 
(Hi) There exists s' S such that Ker{A{s')) — {0}. 

(iv) There exists a constant C ~ C{S) such that \\Ar{s)\\c{Er,Fr) ^ ^ fof 0,11 r > 0. 

Then, if we denote by W{so) the "root subspace" associated to A(sq), that is, the 
linear space generated by the chain of vectors {eo, ei, . . . , e^, . . .} defined as, 

{A-soB)eo = Q, {A- soB)ei= Bea,... [A ^ soB)ek = Beu-i, ■ ■ ■ 

and if we denote by WrisQ^S) the hull of all "root subspaces" associated to Ar{s) for 
all \s ^ sqI < S, s G S , then we have that for S > small enough 

distHiWr{sn,S), W{so)) — > 0, as r — > +oo, 

and therefore there exists a 6 > small such that 

dim{Wr{so, 6)) = dim{W (sq)) , as r ^ +oo. 

Proof See the proof in |ini|3ni[3I]- □ 
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Let us write our operators in such a way that we can obtain the regular conver- 
gence. Consider the following setting. Let E = H'^{R,C'^) and F = ^^(MX^), the 
spaces in the whole real line. Also, = H^{Ir-,'C'^) and Fr = I/^(/r,C^) x C^, the 
spaces in the finite interval Ir and observe that the space F^ has two extra coordinates. 

Define the family of linear operators Pr : E ^ Er and Qr : F ^ Fr as 



and 





V J 



Consider the family of operators A%^,Aq^^, Hoc : E ^ F, defined as, 



+ 



-/ 







+ 



and 



-I 

f'i^oo) -S Aoo 





7rcx>(w)^'cK 



and observe that A^ = ^ + Hoo- Moreover, the operator Aq is a local operator. 
The operator A^ can also be decomposed as 



AS — aO _ „p> 



where 



(5.20) 



(5.21) 



With respect to the operators in a bounded interval, we define, A^jAq^j^t '■ 
Er — >■ Fr as 



Vx 





+ 
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f'i^r) - 
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At- 



0/ 



+ 






u{a{r)) 
\u{b{r))J 



and 



/ \ 



V ) 



and observe that in a similar way, we have 



(5.22) 



with 



Br 



u 




We have the following, 

Proposition 5.8. With the notation above, for any s G {Res > ^v}, we have 
(i) The sequence of operators Aq ^ VQ -converges regularly to Aq ^ as r ^ oo. 
(a) The sequence of operators A^ T^Q -converges regularly to A^ as r ^ oo. 
(Hi) The family of operators A^^, Af. are Fredholm operators of index 0. 

Proof, {i) Let us define the auxiliary operator ^ : iJ^ — >■ which is given by, 



A?, 







( ' 


-A 


Vx 




/'(<i>oo) 







+ 












^ 


0/ 






u{a{r)) 
\u{b{r))J 



where we consider that /'(<i>oo) is restricted to the interval 7^. Notice that Af^^ 



Al^ + Br where 



Br 



f \ 

-./'(*oo) + -Aoo+A 






and observe that \\Br 



f'i^r)\\ 



0, see Lemmas 



4.11 



as r 
and 14.14 



since A^ Aoo and ||/'($oo) 



But from |3] we know that Aq ^ -converges regularly to Aq^. This conver- 
gence is not trivial at all and it uses deep techniques like exponential dichotomy. Also, 
it is worthwile to mention that this result is implicit in the work of Beyn and Lorenz 

m- 

The fact now that ||Sr||£(_Er._F,.) ~^ implies easily the result. 
(ii) Once we have obtained the convergence for the "local" operators and recalling 
that Ar = AQ r + Ilr, A^ = ^0,00 + Hoo wc obtain the VQ regular convergence from 



the VQ regular convergence of Aq j. to Aq ^, the VQ convergence of to He 
the fact that both Hr and Hoo are operators with a 1-dimensional rank. 



and 



(Hi) Let us divide the proof in two parts. 
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A'l is Predholm with index 0. Operator is defined in tlie finite 
interval Ir where wc have the compact embedding H^{Ij.,C) M- Z/^(/r,C). This 
implies in particular that the operator 







( ' 


-A 





















\0) 






0/ 



+ 






u(a(r)) 
\u{b{r))J 



is a compact operator from Ej. to F^., since it is a bounded operator from L'^{Ir,C) x 

L'^{Ir,C) to Fr. 

Hence, the operator : Er ^ Fr is a Predholm operator of index if and only 
if the bounded operator : E^ ^ F^, given by 







is a Fredholm operator of index 0. But this is very easy to show, since Ker{Dr) = 
{{u,v) € Er : u — constant, v = constant} = C x C and therefore dim(Ker{Dr)) = 2. 
Moreover, the rank of Dr is L^{Ir) x L^{Ir) x {0} x {0} C Fr which has codimension 
2. 

(m-2) A^ is Predholm with index 0. Observe that the operator A^ can be 
decomposed as A^ = F^ + K^o + where liao is given as above (and is a compact 
operator since it has rank=l), and the other two operators are given as 



and 







-/ 
V{-)-s 



where the potential V{x) is piecewise constant and it is defined as 

y[X)-<^ /'(I), XG{0,00). 

But the fact that ^od{x) — >■ as x — )• — oo and ^oo{x) — >■ 1 as x — >■ +oo, implies 
that /'($oo(x)) — V{x) — > as a; — > ±oo and therefore, the operator Kr^ '■ Erxi — > -FLo 
is a compact operator. Hence, A"^ is a Fredholm operator of index if and only if 
F^ is a Fredholm operator of index 0. 

The operator F^ is written as 



+ M{x,s) 



with M{x, s) the piecewise constant matrix function, 
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a; < 0, 



and recall that both /'(0),/'(l) < 0. 

To show that is Fredholm with index 0, we will show that Ker{F^) = {0} 
and R{F^) = L^[^,<C'^)- The fact that Ker{F^) = {0} is proved as follows. Let 
{u,v) e iJi(i?,C2) such that 



oo 



U 

vl ^ \0 



Then, if we consider this equation in a; < (reps, a; > 0), it is a linear 2x2 ODE with 
constant coefficient whose solution can be obtained explicitly. Since {u, v) e 77^ (M, C) 
we have that {u,v) is a bounded function as |a;| — oo and therefore, necessarily the 
behavior of the solution as a; — ^ — oo (resp. x — >■ +oo) is completely determined by 
the spectral decomposition of the matrix M^{s) (resp. M+{s)). 

Direct computations show that both matrices M_(s) and M^{s) are hyperbolic 
matrices (no eigenvalues with real part), each of them has one eigenvalue with 
positive real part and the other with negative real part. If we denote by ap{s) the 

eigenvalue with positive real part of M_(s) which has ( "I n ) as its associated eigen- 



vector (unstable manifold of of M_(s)) and by w„(s) the eigenvalue with negative 
real part of A/_|_(s) which has ( ^ ) as its associated eigenvector (stable manifold 



of of A/+(s)) then we necessarily have that 



"SiUc^e-^^f ;,Vfor.<0, f"iiUc„e--f \^ 



for some constants Cp,c„ G C. But since (u^v) G i?^(M,C), we necessarily mut have 



1 \ /I 

"pi 



and this is impossible unless Cp = Cn = since ReQ;p(s) > and Rea;„(s) < and 
therefore both vectors are linearly independent. This shows that {u, v) = (0, 0) and 
therefore, Ker{F^) = {0}. 

To show that R{F^) = L'^{R,C'^), we apply TW, (Lemma 1, p.l37)]. Observe 

that again, the proof of this result uses that both vectors ( ^/ \ ) and ( ^/ \ ) are 



linearly independent and generate the complete space C^. □ 



Remark 8. Behind the proof above we have implicitly used the fact that the 
operator F^ has an "exponential dichotomy" in the luhole real line M. We refer to 
13 \2iA for literature relating exponential dichotomies and Fredholm operators. 



Remark 9. Observe that the proof that Af. is Fredholm of index is valid for 
all s e C, while the proof that is Fredholm of index uses in a decisive way that 
Res > max{/'(0), /'(I)}. This is related to the fact that the essential spectrum of L°° 
is contained in {z £ C : Re z > max{/'(0), /'(I)}}, see Remark^ 
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Theorem 5.9. For every fixed e > 0, there exists an r^ > such that for all 
r > rQ we have a{Lj.) D {Re z > + e} = {s(r)}. Moreover, s{r) < is a simple 
eigenvalue of and s{r) as r +oo. In particular, the unique stationary 



solution of (4.101 is asymptotically stable. 

Proof. Observe first that from Corollary |5.2| we have that there exists i?o > 
large enough and independent of r such that a{Lr) n {Re 2; > —v} C {\z\ < Rq} 
and therefore, the part of the spectrum of with Re 2: > — ;/ is uniformly bounded. 
Hence, from now on in the proof of this theorem, we will only consider s G C with 
\s\ < Rq and Res > — i^. 

Observe that s E <j{Lr) if and only if Ker{A''^) 7^ 0. Moreover, notice that if 
s is such that Res > — i^, then Ker{A''^) ^ {0} if and only if s is an eigenvalue 
of L°° . Hence, from the spectral analysis performed above for we have that 
Ker{A^^) = {0} for all s with Re s > —v except for s = 0, for which Ker{A''^) is one 
dimensional and it is generated by the vector function {•^'aoi'^'Lo)- 

Let us calculate the "root subspace" associated to sq — 0. Following Theorem 
we have that cq — ('f^, To calculate ei, we need to solve A^ei = BooGq, 

lat 



5.7 



where was defined above. That is 



which can be written as 



Ux ~ V — 

Vx + f'i^oo)u + XooV + TT^iu)^'^ = ^'^ 



00 



or equivalently £Jf (u) = u^x + f'{^)u + XooUx = *Jx) - •^oo(m)$L S span($^) which 
has no solution since $^ is an algebraically simple eigenfunction of Hence 

W{Q) = span | | and dim(W^(0)) = 1. 

Therefore, from the VQ regular convergence of A^ to A^^ given by Proposition 
5.8 and applying the definition of Wr{so,5) and the results from Theorem 5.7 we 
have the following: 

i) All values s(r) e C with Res(r) > -1/ and \s{r)\ < Rq for which Ker{Ar'-''^) ^ {0} 
satisfy s(r) — as r — +00. 

ii) There exists a ^ > small such that for r large enough we have dini(Wr(0, S)) = 1. 

In particular s(r) from i) is a real number since if it were a complex number, then 
its complex conjugate s(r) would also satisfy Ker{Ar^^^) ^ {0}, since if A^'r'^^ {u,v) = 
(0, 0) and Im s(r) ^ 0, then Im (m, z;) ^ (0, 0). Moreover, since all coefficients of A^}-^^ 
are real (except for s(r)), we will have A^r^^(u^ v) — (0, 0) and therefore we will have at 
least two numbers, s(r) and s(r) in the set {s : Res > —v, Ker{Af.) ^ {0}}. From 
i) we will have that both of them have to aproach and therefore, dim(Wr(0, 6)) > 2, 
which is a contradiction with ii). This shows that s(r) e M. Hence, s{r) is a real 



eigenvalue of Lr, but from Proposition 5.4 we have that s(r) < 0. □ 



6. Numerical experiments and open problems. In this section we propose 
numerical examples that show the efficiency of the methods analyzed in this paper. 
Notice that the implementation of numerical methods necessarily passes by the trun- 
cation of the interval and the use of certain numerical schemes applied to the truncated 
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equation. The stability analysis carried out for the equation in a bounded domain in 
the previous section shows that if the time is large enough but finite an appropriate 
initial data will be close to the unique stationary solution. Moreover, for this fixed 
time, if the numerical scheme is convergent then for an appropriate refinement of the 
discretization mesh, we will obtain a numerical approximation of the equilibria which, 
in turn, will also be an approximation of the travelling wave. Which is the main goal 
of this paper. 

Observe that at this moment, one step further may be given which consists on 
the analysis of the dynamics of the equations obtained by the discretization, maybe 
discretizing both space and time or just discretizing only the space. The analysis 
of the dynamics of the resulting discrete equations and the comparison with the 
dynamics of the continuous equations is an interesting subject that will be analyzed 
in a forthcoming work. 

We consider the prototypical Nagumo equation 



Ut = + m(1 — u){u — a), a; e M, i > 0, 
u{x, 0) = uq{x), X e M, 

for which an explicit travelling wave solution u{x,t) — $(a; 

1 



a e 



(6.1) 



1 



E/V2' 



V2 



ct) is known, namely 

(6.2) 



X e 



As we have already mentioned in the Introduction, we only consider here monotonic 
increasing waves, taking values at — oo and 1 at +oo, since the other case, that of 
monotonic decreasing fronts, can be dealt with by simply changing variables from x 

to —X. 

It is clear that f{u) = u{\ — u){u — a) fulfils the hypotheses of Theorems 2.1 and 



3.3 and we can use the change of variables (3.1 1 to approximate both the asymptotic 



travelling front $ and its propagation speed. In this way, for fixed J > 0, we consider 



the numerical integration of (1.16) in the interval [—J, J], this is 

m 



v{x,0) = va{x) 
v{~J,t)^0, v{J,t) 



+ v{l - v){v - a), X e[-J,J], t>0, 



LH-J,J) 

X e [-J, J], 

1, t>0. 



(6.3) 



We applied the method of lines to integrate ( 6.3 1 up to time T = 150, for a — 1/4. 
For the spatial discretization, we use standard finite differences formulas, centered for 
the approximation of v^, on the uniform grid Xj = — J + jAx, 1 < j < M — 1, for 
Ace = 2 J/M, and different values of J and M . The nonlocal term is approximated 
by using the scalar product of the vector with the values of Vx at the grid points Xj. 

For the time integration of the spatially semidiscrete problem we use the MAT- 
LAB solver odel5s. Since we are interested in computing both v and A„, it is conve- 



X e [-J, J], i > 0, 



nient to reformulate (6.3) as a Partial Differential Algebraic Equation 

Vt = Vxx + Kvx + v{l - v){v - a), 

= K{vx,vx) +F{1), t>0, 

vix,0) = vq{x), xe[-J, J], 

w(-J,t) = 0, v{J,t) = l, t>0, 
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(6.4) 



The results shown in this section are obtained with the following options for the solver: 
e = ones (M-1 , 1) ; 

Massmatrix = spdiagsC [e;0] ,0,M,M) ; 

options = odesetC Mass Massmatrix, 'MassSingular' , 'yes' ,.. . (6.5) 
'MStateDependence' , 'none' , . . . 

'AbsTol' , [le-8*ones(M-l,l) ; le-10] , 'RelTol' , le-8) ; . 



In the Figures below we show the results obtained for different possible choices of 
the initial data uq{x). These plots show how for J and M large enough the solution 
to the discrete problem seems to converge exponentially fast in time to a stationary 



state close to the stationary state of (1.9). The same happens with the value of the 



propagation speed. These numerical results are in agreement with the ones reported 
in [5l I25j and, up to a great extent, are to be expected from our theoretical analysis. 
However, let us notice that while Theorem |3 . 3| guarantees convergence to the equilib- 
rium of the problem in the whole line for a wide class of initial data. Theorem |4.10| 
guarantees convergence for the problem in a bounded interval only for initial data 
close enough to the equilibrium and in a large enough interval. How "close" and 
"large" are "enough" , is not really specified. 

Several additional questions arise now related with the asymptotic behavior of 



problem (6.4 1 as t ^ oo. For instance, we know that the unique equilibrium of (1.161 



approaches the unique equilibrium of (1.151 as the interval grows to M. A natural 



question now is what the rate of convergence with repect to the length on the interval 
is. With the notation in this section, this is the convergence with respect to J. 

Concerning the boundary conditions, other variants are meaningful and could in 
principle lead to a faster convergence to equilibrium, such us homogeneous boundary 
conditions of Neumann type or even more sophiscated conditions like transparent 



boundary conditions. Notice that then problems (1.9 1 and (1.151 lead to different 
IVBP. 

Another issue is the effect of the numerical approximation. One could address 
for instance the study of the speed of convergence towards the asymptotic profiles 
depending on the chosen numerical scheme, the mesh size, etc. In particular, it would 
be natural to analyze the question of whether upwinding yields better convergence 
rates. The same questions make sense for fully discrete approximation schemes. 

Finally, let us notice how the worst approximation results displayed in Figure |6.4| 



illustrate the importance of capturing properly the front of the asymptotic profile. 



In other words, the importance of controlling the value of the phase xi in (3.8 1. A 
careful study of the dependence of this location on the initial data uq is in order. 

All these questions are beyond the scope of the present paper. 

Example 1. We consider the linear initial data 

a; + J 

Mo(a;) = xe[-J,J]. (6.6) 



The results plotted in Figures [6T] and were obtained with J = 40, Ax — 0.1 and 
Ax = 0.025. 

Example 2. We consider the initial data 

uo{x) = 1 (l + 0.53^ + 0.47 sin ( - |^)) • (6.7) 
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Fig. 6.1. Solution for uq in Example [j| (left) and evolution of A„ (right) for J = 40 and 
Ax = 0.1. 
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Fig. 6.2. Error in the approximation of c for Example^ Left: Ax = 0.1, Right: Ax = 0.025. 
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Fig. 6.3. Solution for uq in Example \^left ) and evolution of Xy (right) for J = 40 and Ax = 0.1. 
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Fig. 6.4. Error in the approximation of c for Example Left: Ax = 0.1, Right: Ax = 0.025. 
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Fig. 6.5. Solution for uq in Example (left) and evolution of (right) for J = 40 and Ax = 0.1. 



Example 3. We consider the initial data 



Uo{x) 



0.2, 
0.8, 



if a; < 
if a; > 0. 



(6.8) 



In this case uq does not satisfy the boundary conditions. We chose this example 
because if we consider the natural extension of uq to the whole line (by 0.8 to the 
right and 0.2 to the left). Theorem 3.3 guarantees the convergence of the solution of 
(pjl) to $ in ([6!2|. 



Appendix A. Technical lemma. We include in this appendix a technical 
result, where we study in detail the behavior at ±oo of the bounded solutions to a 
certain kind of second order differential equations with variable coefficients. 

Lemma A.l. Let c S M and let us consider the second order (non homogeneous) 
scalar differential equation 



ip"{x) + cip'{x) + a{x)ip{x) = f{x), Xq < X < oo. 



(A.l) 



where 
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Fig. 6.6. Error in the approximation of c for Example^ Left: Ax = 0.1, Right: Ax = 0.025. 



i) a{x) is a bounded piecewise continuous function satisfying 



limsup \a{x) - a| < Mie , 



ii) the function f{x) satisfies 



\fix)\<M2e- 



(A.2) 



(A.3) 



where 9,t > 0. Then any bounded solution ij) of (A.l) tends to as x ^ oo. More- 
over, if < oj < min{6, —ri,T}, with ri = -'^-^^^- i^L < Q, then 

\^{x)\, \^'(x)\, \i;"{x)\ < Me-"% for x > X2, 



for some constant M > 0. 

Proof. Let us rewrite the equation |A.1| as 



■0" + c^p' + a?/> = (a - a{x))ij{x) + f{x) := b{x), 



(A.4) 



and observe that from i) and ii) we have \b(x)\ < A'le^^^ for some M > Q and with 
7 = min{r, 0} > 0. 

The roots of the characteristic equation associated to the homogeneous equation 



of (A.4) are precisely 



-c - \/c^ - 4a -c + - 4a 
ri = < and ra = > 0. (A. 5) 



By the variation of constants formula, any solution ^ of (A.4) is of the form 

1 



ri(x-s) _ r-2(x-s) 



)b{s)ds, (A.6) 



for C, D G M. If require further that -0 is bounded, the only possible choice for D is 

1 



D ^ 
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e-''^'h{s) ds, 



(A.7) 



leading to 



But, if ri + 7 7^ 0, then 



ri - r2 



'b{s) ds. 



(A. 



'b{s)ds 



< M I e 



-(ri+7)s ^ 



and if ri + 7 = 0, then 



'b{s)ds 



Xo 



Moreover, since r2 + 7 > 0, we have 



f 

J X 



'b{s)ds 



< M 



,-(r-2+l)x 



Plugging these estimates in ( A.8), and with some simple computations we obtain that, 
if Ti + 7 7^ then 

and if ri + 7 = 0, then 

mx)\ < Cie'-'^ + C2e-^^x - xo) + ^36"^% 
from where the conclusion for tp follows easily. 



To obtain the bounds for 'ip'i^) we just take derivatives in (A.8 1. We obtain an 
extra term, b{x), and the rest of the terms are estimated similarly as in the case of 
il^{x). To estimate tp"{x) we use the equation satisfied by ip and the bounds obtained 
for ^p{x) and ip'{x). □ 
Remark 10. The same conclusions of the previous Lemma hold if we are dealing 
with the interval —00 < x < xi. In this case we need to specify the behavior of the 
functions a(-) and /(•) as x —00 and the conclusion is the exponential decay of the 
solution as X —00. 



Appendix B. Proof of Lemma |5.12[ Finally, we include in this appendix a 



proof of Lemma 5.5 



Proof, (i) and {ii) follow from [191 Section 5.4 and Appendix A]. 

{Hi) Observe first that we know the behavior of $00 (a;), $^(2;) as a; — >■ ±00. 
Notice that the orbit x — )■ (<I>oo(a;), is the heteroclinic orbit connecting (0,0) 

(as X — —00) with (1, 0) (as x — +00) of the ODE, 



U' = V, 

V = -XV - fiU) 



(B.l) 



and therefore the orbit lies in the unstable manifold of (0,0) and the stable manifold 
of (1, 0). Via linearization of the equation in (0, 0) and (1, 0) we can obtain that if we 
define 



ri 



-A- v/A2-4/'(l) 



r2 



-A+ VA2-4/'(0) 



42 



we have 



|$oo(x)|, \^U^)\ < Ce'-\ asx^ -oo, 
|$oo(x) - 1|, \<i>Ux)\ < Ce'-i", asx^+oo 

and therefore 

|/($^)|, I/' ($00 (a;)) - /'(0)| < Ce'■^^ as x ~cx3 

and 

l/(<I'oo)|, Ifi'i'U^)) - /'(1)| < Ce--\ asx^ +00. 
Usmg the equation for $00, that is, = — — f{^oc) = 0, we get also 
< Ce'^^'^, as a; ^ -c5o, \<^'^{x)\ < Ce''^^, as a; -> +00. 



Applying Lemma |A.1| and Remark \10\ we have that if there exists a function w 
then 

w{x)\, \w'{x)\, \w"{x)\ < Ce''2 ^, as X -> -00, 



such that Lg°w = then 



and 

\w{x)\,\w\x)\,\w'\x)\<Ce^ 

where < < r2 and ri < < but arbitrarily close to r2 and ri, respectively. 

Once this estimates have been obtained, we can perform the change of variables 
v{x) = e^^w{x) which will be a function in iJ-^(]R), because of the estimates found 
above for w,w',w". Therefore, v will be a solution of 



(B.2) 



where Xoo{x) — e2^<I)^(x), which is a function in L^(M) because of the exponential 
bounds obtained for ^'^(x) and it is an eigenfunction of the operator v ^ v" + 

+ V associated to the eigenvalue 0. But this operator is selfadjoint 



and therefore there cannot exist a solution of equation (B.2 1. □ 
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